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NOTES ON TOPOLOGICAL INSULATORS 


RALPH M. KAUFMANN, DAN LI, AND BIRGIT WEHEFRITZ-KAUFMANN 

Abstract. This paper is a survey of the Z2-valued invariant of 
topological insulators used in condensed matter physics. The Z- 
valued topological invariant, which was originally called the TKNN 
invariant in physics, has now been fully understood as the first 
Chern number. The Z2 invariant is more mysterious, we will ex¬ 
plain its equivalent descriptions from different points of view and 
provide the relations between them. These invariants provide the 
classification of topological insulators with different symmetries in 
which K-theory plays an important role. Moreover, we establish 
that both invariants are realizations of index theorems which can 
also be understood in terms of condensed matter physics. 
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1 . Introduction 


Topological insulators [41] are new materials observed in nature 
which behave like insulators in the bulk but have conducting edge or 
surface states on the boundary. In physical terms, topological insula¬ 
tors are gapped electronic systems which show topologically protected 
non-trivial phases in the presence of the time reversal Z 2 -symmetry. 
Because of the (odd) time reversal symmetry, topological insulators are 
characterized by a Z 2 -valued invariant. We will use Z 2 to denote the 
group Z/2Z, since this is the notation prevalent in the physics litera¬ 
ture. In the framework of symmetry protected topological (SPT) orders 
[39], the (integer) quantum Hall effect is also viewed as a topological 
insulator in a broad sense, which is protected by a 1/(l)-symmetry and 
characterized by a Z-valued invariant. 

In this paper, we review the topological Z/Z 2 invariants that have 
appeared in 2d and 3d topological insulators and show that they are 
realistic applications of index theory and K-theory in condensed matter 
physics. We stress that the Z 2 index of time reversal invariant topo¬ 
logical insulators can be understood as a mod 2 index theorem. This 
allows us to link together several different representations and thus 
explain them as instances of a common phenomenon. 

For example, the Atiyah-Singer index theorem [7] can be used to 
compute the analytical index of an elliptic operator by its topologi¬ 
cal index. In physics, studying the spectral flow of observables, such 
as counting the net change of quantum charges, can be converted to 
computing topological invariants, e.g. evaluating the Chern character 
in gauge theory. This evaluation can then be seen as an application of 
the index theorem, which of course has intimate relation with K-theory. 

The integer quantum Hall effect gives a hrst example of a topolog¬ 
ical insulator. The hrst Chern number derived from the underlying 
geometry, which used to be called the TKNN invariant [66], describes 
the quantized Hall conductivity and characterizes the integer quantum 
Hall effect. More precisely, the Hall conductance can be calculated by 
the integral of a hrst Chern class over the Brillouin torus T^, which is 
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often compared to the Gauss-Bonnet theorem on surfaces [11]. The 
quantum Hall effect will be discussed from a topological point of view 
in section 2. 

Its cousin, the quantum spin Hall effect [54] was hrst studied in 
graphene, and later observed in HgTe quantum wells with stronger 
spin-orbit coupling [52], More recently, 3d topological insulators were 
also observed in nature [41] . These 2d and 3d topological insulators are 
characterized by a Z 2 -valued invariant, and they are stable under small 
perturbations by impurities or defects, so such materials are called Z 2 
topological insulators. Physically, the Z 2 invariant counts the parity of 
Majorana zero modes, which are quasi-particles slightly different from 
Majorana spinors. 

Since there are many equivalent dehnitions of the Z 2 invariant, one 
main goal of this survey is to show that different descriptions of the Z 2 
invariant can be unihed by a mod 2 index hrst introduced by Atiyah 
and Singer in [9] . The main fact is that taking into account the time re¬ 
versal Z 2 symmetry changes the geometry and topology of the Brillouin 
torus. This is the reason for the appearance of both the topological Z 2 
invariant and the mod 2 index theorem. 

Physicists originally proposed two Z 2 invariants from different per¬ 
spectives. The Kane-Mele invariant was hrst introduced in the study 
of the quantum spin Hall ehect on graphene [44], and generalized to 
3d topological insulators later [33, 55]. The Kane-Mele invariant is 
dehned as the product of signs of Pfafhans over the hxed points of the 
time reversal symmetry. As an alternative to the discrete Pfafhan for¬ 
malism of the Kane-Mele invariant, the Chern-Simons invariant has an 
integral form dehned by the variation of Chern-Simons actions under 
a specihc gauge transformation induced by the time reversal symmetry 
[61]. In the language of index theory, the Kane-Mele invariant is an 
analytical index and the Chern-Simons invariant is a topological index. 
Together they ht perfectly into a mod 2 index theorem. So we view 
them as diherent aspects of the same Z 2 invariant, and we call both of 
them the topological Z 2 invariant from now on. In section 3, we will 
look into the topological Z 2 invariant and its variants from diherent 
points of view. We then show how these dehnitions can be viewed as 
a mod 2 analytical or topological index in section 4. 

In order to model the band structure of a topological insulator, we 
consider the associated Hilbert bundle, usually called the Bloch bundle, 
over the Brillouin torus. When an additional time reversal symmetry 
is present, the Bloch bundle becomes a so-called Quaternionic vector 
bundle, see §3.3.2. This is the reason that Quaternionic K-theory KH 
(or KQ), and with it —via isomorphism theorems— Real K-theory 
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KR, plays an important role in the classification problem of topological 
insulators. More generally, when one considers general finite symmetry 
groups in condensed matter physics, twisted equivariant K-theory can 
be used to classify topological phases with symmetries [31]. A brief 
discussion about the K-theoretic classification of topological phases is 
also included in section 5. We conclude with further discussions of 
an index theorem realizing the bulk-edge correspondence and possible 
generalizations to noncommutative geometry. 

In this survey, we collect ideas and results from both the physical 
and mathematical literature in order to link them —by exhibiting the 
mechanisms under which they are equivalent— and put them into per¬ 
spective under the common aspect of index theorems. We focus on 
those creative ideas and write down the most careful treatment to date 
of those constructions, but at each step provide the necessary setup 
and references where the relevant details can be found. We also in¬ 
clude some original ideas of our own, for example, we derive a formula 
of the Z 2 invariant based on holonomy in §3.2.6. This survey is dif¬ 
ferent from others since we put the topological Z/Z 2 invariants in the 
framework of index theory and K-theory. 

As the literature of topological insulators is already very vast, we 
apologize in advance if we inadvertently missed some material. 


2. First Chern number as the Z invariant 


In this section, we will review the integer quantum Hall effect (QHE), 
the first Chern number, and the bulk-boundary correspondence of QHE 
as an index theorem. The purpose of this section is twofold. First, the 
quantum Hall effect provides the prototype of a topological insulator, 
in which the only topological invariant is given by the first Chern num¬ 
ber. It is helpful to understand the geometric picture of this simplest 
2d model before further generalizations. The second purpose is to fa¬ 
miliarize the reader with the basic language used in condensed matter 
physics. 

In the classical Hall system, a magnetic field B is turned on perpen¬ 
dicular to a plane of finite size. Then one applies an electric current 
and considers a single electron flowing in one direction in the plane, 
say the y-direction. The Hall conductance is computed as the ratio of 
the current density to the electric field. 


where u is the so-called filling factor, e is the elementary charge and h 
is the Planck constant. 
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When the magnetic held is very strong and the temperature is low 
enough, the hlling factor is quantized and takes on positive integer 
values. This phenomenon is called the integer quantum Hall effect, 

( 1 ) (^xy = n—, neN 

The quantum number n, also called the TKNN integer, is successfully 
understood as the hrst Chern number of the Berry curvature over the 
Brillouin torus. 

For an even stronger magnetic held and at extremely low tempera¬ 
ture, the fractional quantum Hall ehect can be observed, in which the 
hlling factor can also take on fractional values. The fractional QHE 
will not be discussed here, since we will only consider Hamiltonians 
with weak interactions in this survey while in the fractional case the 
Coulomb interaction between electrons is dominant. 


2.1. Landau levels. As a warm up, we consider the spectral theory of 
a free Hamiltonian in the Landau gauge in this subsection. Its spectrum 
consists of the Landau levels, which is one of the simplest and most 
important band structures in condensed matter physics. 

Let us hrst forget about the spin property of electrons for the mo¬ 
ment. A free electron in a magnetic held has the Hamiltonian 

Since the magnetic held B is assumed to be perpendicular to the xy- 
plane, we consider the Landau gauge A = (0, Bx, 0) such that B = 
V X A = (0, 0, H), in this gauge 


It is obvious that the commutator [H,py] = 0, so H and Py can be 
diagonalized simultaneously. In other words, in the Hamiltonian py 
may be replaced by its eigenvalue ky, 

^=^[Px + (ky + eBxf] 


Using the so-called cyclotron frequency u = eBim, the Hamiltonian 
is written as 


( 2 ) 


H 


2m 2 mu 


2 
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which takes the same form as the quantum harmonic oscillator. There¬ 
fore, the energy spectrum is given by 



n > 0 


The set of wave functions (eigenspace) with energy corresponding to 
an integer n is called the n-th Landau level. Denoting the harmonic 
oscillator eigenstates by 'ijjn, the wave function of the free electron in 
our model then is expressed as 



(3) 


This state (3) is the first example of an edge state, which by defini¬ 
tion is the eigenfunction satisfying the eigenvalue equation of the edge 
Hamiltonian (2). 

In general, given a bulk Hamiltonian y) which has transla¬ 

tional invariance along one direction —say the y-direction— an edge 
Hamiltonian is obtained from the bulk Hamiltonian by a partial Fourier 
transform. The momenta ky in y direction can then be treated as a 
parameter, so that the edge Hamiltonian becomes a one-parameter 
family ky) as do the energy levels and the eigenfunctions, he. we 

have the following family of equations: 


ky)(l)j^(^X^ ^y) ^niky)4^n^X^ 


parameterized by ky. 

In topological band theory, the band structure, i.e., the Landau lev¬ 
els, defines a complex vector bundle over the momentum space, more 
details will be given in the next subsection. In the simplest case, the 
band structure of a free Hamiltonian gives rise to a trivial vector bun¬ 
dle. 

2.2. Bloch bundle. Thouless, Kohmoto, Nightingale and den Nijs 
(TKNN) [66] first showed that the quantum number n in the Hall con¬ 
ductance (Jxy = ne^/h, which is called the TKNN invariant by physi¬ 
cists, can be topologically identified with the first Chern number. In 
this subsection, we will review the first Chern number of the quantum 
Hall effect based on topological band theory. 

For the 2d square lattice F ~ C its Pontryagin dual is the 
2-torus T^, and it will be called the Brillouin torus from now on. In 
condensed matter physics, a space dual to the original lattice is called 
the Brillouin zone B. In physics parlance, this is the momentum space 
in which the crystal wave vector k lives. 


NOTES ON TOPOLOGICAL INSULATORS 


7 


For electrons in a crystal, the electronic states are described by Bloch 
wave functions 

(4) u{r + R) = M(r), V R G F 

In the following, we use the momentum representation. The quasi- 
periodic eigenstates Mn(k) are by dehnition solutions to the eigenvalue 
equation, 

R'(k)M„(k) = En{'k)un{'k) 

where En{k) is called the n-th band function. We assume that the 
£'„(k) are non-degenerate everywhere (i.e., no band crossing). In gen¬ 
eral, the energy spectrum {F^„(k),n > 1} is referred to as the band 
structure of the model. 

Mathematically, the Mn(k) are sections of a Hilbert bundle E ^ B 
whose hber over k G R is given by the Mn(k) with hxed k. Here one 
uses sections and the compact-open topology; see e.g. [31] for a 
detailed study. All the sections form a Hilbert space, which will be 
called "H. The sections Mn(k) with hxed n form the sub-Hilbert space 
T-in and V, = 

In the physical literature sections are called physical states and are 
written as |M„(k)), or |n,k) to emphasize the quantum numbers. Fur¬ 
thermore they are always assumed to be normalized so that (n, k|n, k) = 

1 . 

There are two types of bands: the occupied bands, literally occupied 
by the electrons, and the empty bands. The valence band is the top 
one among the occupied bands and the conduction band is the bottom 
one among the empty bands. The band gap is dehned as the energy 
gap between the valence and conduction bands. The Fermi energy Ep 
is dehned such that in the ground state all single-particle levels with 
smaller energy are occupied and there is no occupied level above this 
energy. Loosely, it is the top of the collection of occupied electron 
energy levels. It coincides with the chemical potential at absolute zero 
temperature [3] , which is the most precise dehnition. What the position 
of this energy inside the band structure is, depends on the specihc 
material. For insulators Ep lies in a band gap, which we will now 
assume. 

The construction of the bundles is valid for generalized Brillouin 
zones such as or T'^, but in the rest of this subsection we only use 
for concreteness, which is the Brillouin zone in the 2d square lattice 
quantum Hall model. Let us consider the rank one (aka. line) bundle 
for the n-th occupied band over the Brillouin torus. 

(5) T^, 
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This is the eigenbundle to the n-th energy band En, which is assumed 
to be below the Fermi energy. Notice that since the eigenvalues of the 
Hamiltonians are real, there is an order on them. Its sections are the 
n-th eigenstates |n„(k)) : —> £n;k h-)■ nn(k). Assume that there 

exist N occupied bands, then the Bloch bundle is defined to be the 
rank N bundle given by the Whitney sum, 

B = ©li-Cn ^ 

The Bloch bundle is a finite rank Hilbert bundle, whose space of 
sections is the physical Hilbert space. The sections themselves are 
Bloch wave functions. 

The Berry connection is then dehned by a covariant derivative on 
the Bloch bundle [31], whose local form on the component is given 
by 

(6) a"'(k) = i(n, k|(i|n, k) 


so the Berry curvature on has a local expression 
(7) f'^{k) = i d{n,k\ A d\n,k) 

The integral of the Berry curvature over the Brillouin torus gives the 
first Chern number. 


( 8 ) 


r = / dk^ 


day 

dk^ dk,. 


For the Bloch bundle H —)■ T^, the total hrst Chern number is then the 
sum 


N 


Cl = 




n=l 


The generalized Bloch-Hamiltonian for the quantum Hall effect (QHE) 
in a magnetic held in the Landau gauge given by (0,ei?a;) is [66]: 

(9) 

H{ki,k2) ^ + Mi)^ + Ei-ik-A + tik2 -eBxf + U{x,y) 

2m dx 2m dy 

with the potential 


U (x, y) = Ui cos(27rx/a) -|- U 2 cos{27iy/b) 

and boundary conditions Uk^k 2 {x + qa, y) exp(—27ripi//6) = Uk^k 2 {.x, y -1- 
b) = cit rational hux p/q, see [66] for details. In particular the 

action of is restricted to a bigger sub-lattice, such that the non¬ 
commutativity is lifted. As was shown in [66] and in the following 
interpretations [14, 16] the quantization of the Hall conductance can 
be explained as the Chern class of a line bundle as discussed above. 
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Fixing such a line bundle, for instance with the hrst Chern number 
n = 1 is frequently called a quantum Hall state. 

If this type of state is put next to an insulator, which has hrst Chern 
number n = 0 then at the edge there are localized chiral edge modes 
[41], This produces so called helical edge states as these states propa¬ 
gate towards the edge, where they get rehected, since there is a different 
topology in the neighboring bulks, but cannot penetrate deep into their 
own bulk, again because of the topological obstruction, see Figure 2 of 

[41]. 

2.3. Dirac Hamiltonian. The quantum Hall system can be viewed 
as a topological insulator, whose effective Hamiltonian is given by the 
Dirac Hamiltonian. A phase transition is then controlled by some pa¬ 
rameters in the Dirac Hamiltonian. In order to illustrate such phenom¬ 
enon, we analyze 2 -band models generalizing the monopole (see below) 
that depend on a function h : H —)■ where B is usually the Bril- 

louin zone and in our case, to be concrete, the Brillouin torus T^. An 
element k = (fcj., ky) G is called a crystal wave vector. A particular 
choice of h, which we discuss, exhibits a phase transition and will also 
appear again in the study of the quantum spin Hall effect. 

Consider the following Hamiltonian depending on k G 

(10) H = kcr 

where the Uj are the Pauli matrices. This Hamiltonian is non-degenerate 
outside k = 0. Thus restricting it to C we are in the situation of 
the previous subsection, with the Bloch bundle being the sum of two 
line bundles Li (B L 2 . It is a nice calculation to show that the Chern 
numbers are ±1 in this case [14, 18]. The degenerate point at k = 0 is 
thought of as a monopole which is responsible for the “charges” given 
by the hrst Chern numbers. 

There is a generalization to a A^-band model, where now the at are 
spin S = matrices. In this case the Chern numbers are n = 2m 
with m = —S', —S -|- 1,..., S' [14, 18]. 

More generally, consider the 2-Ievel Hamiltonian 

(11) H = hx{k)ax + hy{k)ay + hz{k)az = h(k) ■ cr 

Here cr = [ax, Uy, Uz), where the ai are the Pauli matrices. 

Assume that H is nowhere degenerate, i.e., hi[k) 7 ^ 0, the hrst Chern 
number in this case is a winding number. This goes back to [18] and is 
explained in detail in [25] Appendix A. To see this, dehne the function: 

h(k) 

|h(k)| 


d : ^ d(k) 
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The winding nnmber of d compntes the first Chern nnmber of the 
lowest band, aka. the ground state. 

(12) ci = ^ f {dk,d X dfc d) ■ d dk 

dvr J ^2 

This is the pull back via d of the 2-sphere winding number around 
the monopole of if = k - cr at k = (0, 0, 0), see [18, 65]. This also works 
for general Brillouin zones. The hrst Chern number of the other band 
is then — ci since the sum of the two bands is a trivial rank 2 bundle. 
The particular choice of h^k^, ky) yields: 

(13) 

H(k) = Asm{kx)<Jx + Asm{ky)ay + [M — B{sin^{kx/2) +sm‘^{ky/2))]az 

where A, B, M are parameters. M can be viewed as a mass parameter, 
and B plays the role of a magnetic held. This Hamiltonian appears 
in the analysis of an electron on a 2d square lattice with spin-orbit 
interaction. 

It has the hrst Chern number of the lowest band determined by the 
signs of parameters [64], 

(14) Cl = - ]^[sgn{M) + sgn{B)] 

Using the usual hrst order approximation as k —)■ 0, and setting 
B' = |i?, the Hamiltonian takes the form: 


(15) H{k) = Ak^a^ Akyay + [M - B'{kl + /c^)]cr^ 

which is commonly called a Dirac Hamiltonian. Note that comparing 
this expression to a Dirac operator, the above Hamiltonian has a qua¬ 
dratic correction term B'k^a^- It also appears in the analysis of the 
quantum spin Hall ehect in the quantum well, see §3.1.3. 

This system depends heavily on the relation between M and B, since 
this determines the coupling to cr^. he. M and B determine the locus of 
M — Bk? = 0. If we change from B > M to B < M, then the locus will 
change accordingly, which gives rise to a phase transition [64]. Using 
this point of view, the quantum Hall system manifests itself as the 
simplest topological insulator, as we will now discuss. 

2.4. Bulk-edge correspondence. We have seen that the hrst Chern 
number characterizes a bulk state in the quantum Hall effect. When 
considering an interface with an insulator, an edge state appears. In 
order to better understand this phenomenon, we consider two copies 
of the model (15), but with different ground state Chern numbers, put 
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next to each other along an interface. This is an example of a topo¬ 
logical insulator with a domain wall. We discuss how the connection 
between edge and bulk states can be understood in terms of an index 
theorem. 

For this we simplify the situation to consider a two-band system. 
This is enough to capture the phenomenon as we only need to con¬ 
sider two bands separated by a gap in which the Fermi level lies. The 
Hamiltonians are then of Dirac Hamiltonian type and are generalized 
to include terms that may depend on higher derivatives of x and y, see 
below for an example and [36] for details. 

The generalized index theorem for this type of Hamiltonian connects 
the analytical index (aka. the spectral flow) of the edge states and 
the topological index (aka. the hrst Chern number) of the bulk. The 
spectral flow of the edge states and the variation of hrst Chern numbers 
specifying the bulks establishes a so-called bulk-edge correspondence. 

We illustrate this for the Dirac Hamiltonian (15), whose position 
representation is: 

(16) Hb = -id^cCTx - idyay + {M + Bdl + Bdl)az 

where the subscript b stands for “bulk Hamiltonian”. As mentioned 
previously, the corresponding edge Hamiltonian is obtained by a partial 
Fourier transformation in the ^/-direction; see end of §2.1. 

(17) He = -idxCTx kyay + {M + Bdl “ Bkl)a^ 

One introduces a domain wall at a; = 0 separating two bulk sys¬ 
tems, by letting the mass term become dependent on x. This M{x) 
is assumed to have stabilized limits, hm 2 ._,.±oo ^( 2 ^) = The edge 
eigenstates are assumed to satisfy the eigenvalue equation of the edge 
Hamiltonian 

77g(x, ky^(j)ji(^xj ky^ £'^(/cy)(^ 72 (a:, ky^ 

We make the usual assumption, that as \ky \ gets big, then eventually 
£n{ky) will lie in the conduction or valence band. An edge eigenstate 
is called an edge state, if it connects the valence band and the 
conduction band. This means that as one varies ky, en{ky) hits both 
the valence and the conduction band. Hence for any edge state (f>n{ky), 
Sniky) intersects with the Fermi level Ep = 0 (after a possible shift). 

We let sf{(j)n) be the intersection number of Sniky) with Ep = 0, i.e. 
the spectral flow of the edge state (f)n- If the intersection is transversal 
and there are only hnitely many points of intersection, i.e. points where 
Sniky) = 0, then the spectral flow of the edge state 0^ is 


(18) 


+ -n_ e {-1,1} 


■s/(0n) = n. 
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where n+ is the number of times that the sign of Sn changes from 
negative to positive and n_ is the number of times that the sign of Sn 
changes from positive to negative when increasing ky. Note that sf 
makes sense for all the edge eigenstates, but is 0 for the ones that are 
not edge states. 

Set 

z/ = ^ s/(0n) = N+-N_ 

n 

where the sum is either over all eigenstates or over the edge states, 
Nj^ is the number of edge states going from the valence to the con¬ 
duction band as ky increases and N_ being the number of states going 
from the conduction to the valence band. In order to identify the spec¬ 
tral flow with an analytical index, one introduces an extended Dirac 
Hamiltonian [36] 

(19) He = -iVTidky + T2He 

where Xj are the Pauli matrices and n is a constant parameter repre¬ 
senting a velocity. This is a block 4x4 matrix Hamiltonian with the 
diagonal blocks being zero. The anti-diagonal blocks represent the two 
processes corresponding to 7V+ and iV_. In this form it looks like the 
Dirac operator in the Atiyah-Patodi-Singer index theorem [6]. Based 
on the chiral symmetry of He {{He-, Xa} = 0) one can compute that the 
analytical index of He equals the negative of the spectral flow of the 
edge states, 

inda{He) = -V 

The topological index of He was also computed explicitly [36] by 
investigating the divergence of a chiral current of He, see loc. cit. for 
details. The topological index turns out to be the difference of the first 
Chern numbers specifying the two bulk ground states separated by the 
domain wall, 

-indt{He) = Ci(M+) - Ci(M_) 

By the Atiyah-Patodi-Singer index theorem, the analytical index 
equals the topological index for He and hence by the generalized index 
theorem 

(20) 1 /= Ci(M+) — ci(M_) 

which can be interpreted as the bulk-edge correspondence between the 
two bulk states and the edge states. A physical picture is as follows: 
by stacking two systems with topologically distinct ground states, an 
edge state is created at the interface. The transport property of the 
edge state is then determined by the topology of the junction. This 
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would also work for edge states in the quantum spin Hall effect, see 
e.g. [41], 


3. Z 2 TOPOLOGICAL INSULATORS 

In this section, we review the different dehnitions of the Z 2 invariant 
for time reversal invariant free (or weak-interacting) fermionic topolog¬ 
ical insulators. 

For 2d and 3d time reversal invariant topological quantum systems, 
Z 2 invariants characterizing topological insulators were introduced in 
[44, 32, 33, 55]. We will hrst review their dehnitions in this section. 
There are actually several natural dehnitions for these, using Pfafh- 
ans, polarizations, determinant line bundle and holonomy etc. Besides 
giving these dehnitions, we explain how they are equivalent. 

Physically, the Z 2 invariant can be interpreted as the parity of so- 
called Majorana zero modes. In particular calling the invariant z/, for 
a 2d quantum spin Hall system, u = 0 (resp. 1 / = 1) when there exist 
even (resp. odd) pairs of helical edge states. Note that one helical 
edge state consists of two chiral edge states. Similarly for 3d time 
reversal invariant topological insulators, u = 0 (resp. u = 1) when 
there exist even (resp. odd ) numbers of Dirac cones, i.e., conical 
singularities produced by surface states at the hxed points of time 
reversal symmetry. 

We will go through the details in several subsections. §3.1 is an intro¬ 
duction to time reversal symmetry and Z 2 topological insulators. We 
provide examples of two 2d and one 3d topological insulators. We then 
go on to dehne the 2d Kane-Mele invariant in its original version using 
Pfafhans and determinants, being careful to dehne the setup rigorously. 
We also discuss its 3d generalization in §3.2. Here we hrst consider the 
general structure of Bloch bundles with time reversal Z 2 symmetry. 
The next subsection §3.3 then gives a geometric construction of the Z 2 
invariant using bundle theoretic means. While §3.4 reviews the homo- 
topy theoretic characterization of the Z 2 invariant that leads over to 
the considerations of §4. 

3.1. Background and examples. We hrst recall the time reversal 
symmetry, which introduces Kramers degeneracy and changes the ef¬ 
fective topology and geometry of the game. Next, we review the ehec- 
tive Hamiltonians of Z 2 topological insulators, which can be used to 
dehne elements in K-homology in the future. 

3.1.1. Time reversal symmetry. A physical system has time reversal 
symmetry if it is invariant under the time reversal transformation T : 
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t HA —t. In quantum systems the time reversal symmetry is represented 
by a time reversal operator 0. By general theory [71], 0 is necessarily 
anti-unitary and can be realized as the product of a unitary operator 
and the complex conjugation operator. Again by general theory, in 
momentum space (or the Brillouin zone) T acts as k ha —k. 

A time reversal invariant model is required to have [H{r), 0] = 0, or 
in the momentum representation 

(21) QH{k)Q-^ = H{-k) 

For example, in a two-band system the time reversal operator 0 is 
dehned by 

( 22 ) 0 = iayK 

where ay is the imaginary Pauli matrix and K is the complex conju¬ 
gation. 0 is indeed an anti-unitary operator, i.e., for physical states V’ 
and 0, 

(23) (00, 00) = (0,0), 0(a0-|-60) = a00-1-600, a, 6 G C 

For a spin-^ particle, such as an electron, it also has the property: 

(24) 02 = -1 

This results in the so-called Kramers degeneracy, which is the fact 
that all energy levels are doubly degenerate in a time reversal invariant 
electronic system with an odd number of electrons. In fact, 0 and 0(0) 
are orthogonal: 

(25) (00, 0) = (00, 000) = -(00, 0) = 0 
In general, 0 is skew-symmetric in the sense that 

(26) (00, 0) = (00, 000) = -(00,0) 

The Dirac Hamiltonian H{k) = k^a^ + kyay + {M — Bk‘^)az is 
not time reversal invariant. In contrast, the models for quantum spin 
Hall effect with spin-orbit coupling given in the next section are time 
reversal invariant systems. 

3.1.2. Effective Brillouin zone. The association k —)■ H{k) is formally 
a map H : T'^ ^ Hernin, where HerrUn is the space of Hermitian 
n X n matrices. If we take into account the action of T on T'^, then to 
specify a map with time reversal symmetry it is enough to know the 
restriction of H onto a fundamental domain of the T action, by (21). 
Such a choice of fundamental domain will be called Effective Brillouin 
Zone (EBZ). In the 2d-case such an effective Brillouin zone is given by 
a cylinder C = x I G and in 3d by x / C T^, see Figure 1. 
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Figure 1. Effective Brillouin zone in 2d and 3d 


If one considers the boundary of this EBZ, then it still has a non¬ 
trivial action of T and we can consider a fundamental domain of the 
boundary of the EBZ, which we call . In physics parlance, this 
region characterizes the edge states using a unique coordinate. For 
example, = [0,7r] in the 2d effective Brillouin zone [0,7r] x [—vr, vr] 
and is the square [0, tt] x [0,7r] in the 3d case, see Figure 1. 

Although the usual Brillouin zone is T'^, the spheres play a special 
role. First they can be viewed as local invariants as in Simon’s interpre¬ 
tation of the TKNN integers [65] or generally as enclosing monopoles or 
higher codimension degeneracy loci. For instance, from [68], we know 
that for a family of Dirac Hamiltonians i/ G i5, its eigenvalue degen¬ 
eracy is generically of codimension three, so in 3d for a generic family 
each point k corresponding to a degenerate eigenvalue, aka. monopole, 
is enclosed by an 5^. Secondly, the spheres occur through effec¬ 
tive maps as in the Moore-Balents description below §3.4, or the more 
complicated arguments presented in §5.2. Finally, they appear in the 
stable splitting, see §5.4. In particular, they arise through pull-back 
via the collapse map g : T'’* —)■ S''^, see loc. cit.. 

Thus, we will look at maps H : X ^ HerrUn where X can be T'^, S''^, 
or Depending on the situation, the target may be a subspace of 

all Hermitian matrices, see §3.4. 
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3.1.3. Quantum spin Hall effect. We now turn to the Hamiltonians of 
the quantum spin Hall effect, which appears in both graphene and the 
HgTe quantum well. These are 2d time reversal invariant models. Hal¬ 
dane [40] first investigated the quantum Hall effect in graphene which 
has a honeycomb lattice structure. The Kane-Mele model [45] fur¬ 
ther took spin-orbit coupling into account and generalized the spinless 
Haldane model. They hrst predicted the quantum spin Hall effect in 
graphene with the effective Dirac Hamiltonian 

5 5 

(27) Hkm = dt(k)r; 4- 

1=1 l=i<fj 

where for i = 1,... ,5 the gamma matrices Tj are dehned via 

r = (da, 0 So, cr^ 0 So, (Ty 0 Sx, CTy 0 Sy, ay 0 sfj 

where cxj are as before and the Sj are Pauli matrices as well, representing 
the electron spin. Their commutators dehne the T^: 

The time reversal operator in this model is defined by 


(28) 0 = i(cTo 0 Sy)iP 

It is straightforward to check that 

0r,0-i = Ti, 0r,,0-^ = -Ti,- 
Combining this with the properties of the coefficients [45] , 
di{-k) = di{k), dij{-k) = -dij{k) 


one obtains that the Hamiltonian Hkm is time reversal invariant. 

The coefficients also depend on a parameter, which is called the 
nearest neighbor Rashba coupling constant. If the Rashba coupling 
constant vanishes, then the Hamiltonian above is invariant under spin 
rotations along the z-axis, i.e., [Hkm^Sz] = 0. Then it splits into two 
components for the up and down spins. 


H 


t«t(k) 0 \ 

I 0 


where 


(29) 


i7-f- — diax — di2<Ty + (<^2 — di5)az 
Hi = dia^ — di2<Ty -l- (<^2 + dif)az 


This system is basically comprised of two copies of the Haldane model. 
Each component from the above dehnes a quantum Hall system, which 
is characterized by its first Chern number. The total Chern number of 
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the Kane-Mele model is always zero, since for a time reversal invariant 
ground state ci^ = —ci; and hence 


Cl = cit + = 0 


It is, however, possible to dehne a new topological invariant, i.e., the 
spin Chern number 


(30) 



In general, the Rashba coupling is non-zero, then the spin-up and 
spin-down components are coupled together and the dehnition of the 
spin Chern number is no longer valid. However, there always exists a Z2 
invariant for the quantum spin Hall effect, for example, the Kane-Mele 
invariant discussed in section §4.1. 

If a topological insulator can be specihed by a Z2 invariant, we call 
it a Z2 topological insulator. 

The quantum spin Hall effect was also predicted by Bernevig, Hughes 
and Zhang [17] for the HgTe/CdTe quantum well with stronger spin- 
orbit interaction, which was observed by Konig et ah [52] soon after 
the prediction. The effective Hamiltonian for this so-called quantum 
well is given by 



where 


(31) H(k) = Akx(Jx + Akyay -|- (M — -|- (C — Zlk^)(To 


here A, B,C,D are expansion parameters and M is a mass or gap 
parameter. Notice that this effective Hamiltonian, up to the crQ-term 
—which actually drops out in the quantum Hall response—, is of the 
form discussed in §2.3. This system again decouples, so that the simpler 
version of the Z2 invariant given by Cspin is applicable. 

3.1.4. 3d topological insulators. Later on, 3d topological insulators were 
also observed in nature. In 3d the role of edge states is played by 
codimension-1 surface states. The observed surface states look like 
double cones in the spectrum or energy in the momentum dispersion 
relation in physical terms. They take the form of local conical singular¬ 
ities + y'^ = and are called Dirac cones by physicists, see e.g. [46]. 
A point in the Brillouin torus over which a Dirac cone lies is called a 
Dirac point. In 3d topological insulators, a Dirac cone is sometimes 
also referred to as a Majorana zero mode. Note that a Majorana zero 
mode should not be confused with a Majorana fermion. 

If the system is invariant under time reversal symmetry then the 
so-called Kramers degeneracy forces level sticking at any hxed point of 
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the time reversal symmetry, as discussed in the next section. If these 
singularities are not extended and are conical one has Dirac cones. The 
number of Dirac cones determines whether the system is a topological 
insulator, that is, if a material has odd number of Dirac cones, then it 
is a stable 3d topological insulator. 

For example, the effective Dirac Hamiltonian for Bi 2 Se^ with spin 
orbit coupling has a single Dirac cone [74]. The Hamiltonian is 

H{k) = Aikxax^Sx + Aikyax®Sy + A2kzO'x®Sz + M{k)az^So + eo{k.)l4 

where /4 is the 4x4 identity matrix, 

, M(k) = M- B,kl - B2(kl + kl) 

' ' eo(k) = C + Dikl + D^ikl + fcj) 

and the Ai,Bi, C, Di and M are model parameters. 

3.2. Kane—Mele invariant. In this section, we carefully present the 
relevant structures that the time reversal symmetry induces on a Bloch 
bundle and then use these structures to dehne the Kane-Mele invari¬ 
ant in 2d and 3d, following their original construction and hnish with 
a holonomy interpretation of the invariant. We will also discuss a co¬ 
ordinate free bundle theoretic description of these constructions based 
on canonical orientations in the following section. 

3.2.1. Kramers degeneracy. In this subsection, we discuss the geomet¬ 
ric setup for time reversal invariant topological insulators. We are fol¬ 
lowing physical notations, this means that we use sections. These are 
allowed to be non-zero only locally, or in other words locally dehned. 
So we have to stress that in order to get an global object, sometimes 
extra work is needed to glue local sections together. 

Consider a time reversal invariant Dirac Hamiltonian 77(k), k G 
T'^, whose eigenvalue equation is i7(k)M„(k) = i7„(k)n„(k). For each 
physical state |M„(k)) = |n,k), 0|n,k) gives another state with the 
same energy, and such doublet is called a Kramers pair. At a time 
reversal invariant point k = T(k) = —k G T'^ time reversal symmetry 
forces degeneracy or a level crossing. That is both |M„(k)) and 0|Mri,(k)) 
are states with momentum k = T(k) = —k and the same energy. 
They are linearly independent, due to (25). This is called Kramers 
degeneracy. 

Starting with a lattice model, in the reciprocal lattice the time re¬ 
versal invariant points are given by 

= . K, 


—Fj = Fj-f iij ■ K, i.e. Fj 
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where K is the reciprocal lattice vector and the coordinates of rii are 
0 or 1. For example, in a 2d sqnare lattice with K = (27r, 27r), we have 
fonr time reversal invariant points F = {(0, 0), (0, tt), (tt, 0), (tt, tt)}. So, 
for instance, the Brillonin torns has fonr hxed points and similarly 
has eight fixed points of the time reversal symmetry. 

Hence in snch a system, the Bloch bnndle does not split into line 
bnndles according to the energy eigenvalnes, bnt nnder the cnstomary 
assnmption that the Kramers degeneracies are the only degeneracies, 
the Bloch bnndle does split into a direct snm of rank 2 vector bnndles, 
as we now discnss. We assnme from now on that indeed all degeneracies 
that are present are forced by Kramers degeneracy. Physically, this can 
for instance be dne to a spin-orbit conpling. First, let’s assnme that 
the Bloch bnndle has even rank 2N, so that there are 2N states, which 
at the hxed points form pairs. Near a hxed point, the pairs may, and 
nnder the assnmption will, split. So, in a neighborhood of a hxed point, 
each Kramers pair will give rise to two energy levels, spanned by say, 
|w 2 n-i(k ))5 |w 2 n(k)). We will nse the notation: 

(33) |Mn(k)) = |M2n-l(k)), K(k)) = |M2n(k)) 

In this case, if there are only hnitely many occnpied states, there will 
be no crossings between levels from diherent pairs (33).^ This implies 
that if the Bloch bnndle has even rank 2N it can be split nsing the 
Kramers donblets to dehne rank two snb-bnndles 

s.t. H = 

By the same type of argnment if the Bloch bnndle has rank 2N+1, 
then it splits into rank two bnndles and one line bnndle. We will not 
consider this case farther as one can split oh the line bnndle from the 
Bloch bnndle. 

The bnndle B and the snb-bnndles also carry the operation of 
0 which is compatible with the operation T on the base. As we will 
discnss in §3.3.2 this extra strnctnre makes B and the Vn into so- 
called Qnaternionic vector bnndles. Sometimes, we will wish to em¬ 
phasize this and then we will write (H, 0) = 0)„. Notice 

that 0 indnces an isomorphism H ~ H. More precisely we have that 

r*0(v„) ~ K. 

quick proof goes as follows. Suppose that for some n the levels 2n and 2n-|- I 
cross. Then this must happen at a time reversal invariant point. This implies that 
at this point the levels 2n and 2n -I- 1 cross for all n. This will leave the first levels 
and last levels unpaired. Notice that in an infinite system, say indexed by Z, such a 
cross over from the pairing 2n—1, 2n to 2n, 2n+l can happen. This observation was 
the initial spark for the definition of the Kane-Mele Z 2 invariant and is pertinent 
to the discussion of the edge state and the bulk-boundary correspondence. 
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One consequence of this is that all the odd Chern classes of (S, 0) 
are two-torsion 2c2i+i{B) = 0 G Z), since B B implies 

that Ci{B) = CiiB) = (—l)*Ci(f3). Hence they vanish when extending 
the coefficients to M as is done under the deRham isomorphism or if 
the cohomology is torsion free. Real coefficients are always used in 
the computation of Chern classes via integrals, i.e. Chern-Weil theory. 
In particular this shows that the Ci{B) = 0 G for instance 

when computed using integrals over curvature or holonomy. Thus the 
first non-trivial information will be contained in the second Chern class 
C 2 {B) G or more generally the of the base of the bundle. In¬ 

deed, for 4-dimensionaI time reversal invariant models the fundamental 
topological invariant is the second Chern number [61]. However, 
is possibly non-zero only for manifolds of dimension greater or equal 
to 4, so it is not useful directly in 2d and 3d. This is the reason for 
the introduction of the Z2 invariant. For 3d one can use dimensional 
reduction from 4d to 3d to understand how the Z2 invariant appears 
as the dimensional reduction of C2. 

The vanishing of ci{B), however, also implies that the first Chern 
class of its determinant line bundle vanishes, i.e., ci{DetB) = 0 G 
since the cohomology of is torsion free. In addition, this 
determinant line bundle has a global section, which is one of the main 
ingredients of the Z2 invariant. The main key to understanding the Z2 
invariant geometrically can be traced back to [32] , which we now make 
precise. 


3.2.2. Geometric setup. One rough idea about how to obtain the in¬ 
variant is given as follows. At any fixed point ko, ko = —ko, 0 is a 
skew-symmetric (26) anti-linear morphism on the fiber over k. After 
choosing a basis, a natural object to study is the Pfaffian of the result¬ 
ing skew-symmetric matrix. Since we can reduce to rank 2 bundles 
this is just given by the top right corner of a 2 x 2 matrix. The actual 
choice of basis plays a role, as the sign of Pfaffian changes if for in¬ 
stance the two basis vectors are permuted. The Kane-Mele invariant is 
a way to normalize and compare these choices of Pfaffians at all fixed 
points simultaneously. For this one needs to make a coherent choice of 
normalization and show that the result is independent of the choice, 
which is what the invariant does. A key fact that is used is that the 
determinant line bundle has global sections. This can be seen in several 
different ways as we now discuss. 

To study the choice of basis, first observe that at any fixed point ko, 
ko = — ko, 0|0) and |0) are linearly independent by (25), so that we 
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can choose 0|M^(ko)) and 0|M^(ko)), s.t. 

0|Wn(ko)) = \ul(ko)) and 0|M^(ko)) = -|M^(ko)) 

where the second equation follows from the first by 0^ = —1. 

In contrast if |M^(k)) and |M^(k)) are states, i.e. local sections, defined 
on a subset V that is stable under T, i.e. TiV) = V, and contains no 
fixed points, i.e. VnFix(T) = 0, then, since B is trivial over T‘^\Fix(T) 
by the assumption that there are only degeneracies at the fixed points, 
and 0 commutes with hf, we get that away from a fixed point 

(34) u;(-k) = ul(-k) = e'>^«<-'‘>euS(k) 

the relation 0^ = — 1 implies that e*^"’ 

Together these two facts imply that the do not extend as contin¬ 
uous sections to the hxed points in any T stable V that does contain 
fixed points. However, one can cobble together a section using both 
and on certain 1-dimensional domains, for instance in the case of a 
1-dimensional Brillouin zone. 

In particular, let H be a small 0 stable path through a hxed point 
ko. Decompose V as H V~ with T{V^) = V~ and fl V~ = ko. 
Now we can £x sections |M,'j(k)) and |M^(k)) independently on V~^ and 
V~, by continuing the states continuously to the fixed point ko from 
both sides. By perturbation theory, we can choose these to match up 
if necessary after multiplying them with a local phase factor. 

lim M^(k) = lini ul{-k) and lim ^^(k) = lini M^(-k) 

k—^ko —k^ko k—^ko —k—^-ko 

The reason for this is that after subtracting a constant term, the 
family is given by hf(k) = x(k) ■ cr, with x(ko) = 0 and x(k) an 
odd function of k. This means that up to rotation, linearly near kg: 
H(k) = (k — ko)a^ -|- O(k^), whence the claim follows by standard 
perturbation theory. Here we need that we assumed that the fixed 
points are isolated. 

These equations mean that one can glue two sections together over V, 
which we call Iw^fk)) and Iw^^fk)) with |MDk))|y+ = ItiDk)), |MDk))|y- 
\ul{k)) and |M"(k))|y+ = \ul{k)),\ui\k))\v- = K(k)). Likewise, in 
this case, we can glue together the functions x^(—k) and x^(k) at ko 
to form one function Xn{k). For these sections, we obtain that 

(35) ul^{—k) = and m^(— k) = — 

Strictly speaking, to follow the literature, we only need the 1-d ver¬ 
sion for the base with involution T(k) = —k with k G [0,27r] 
mod 271. In this case we have fixed points 0,7r and V~^ = [0,7r], V~ = 
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Figure 2. Sections in the Brillouin zone 


[—71,0]. We fix continuous normalized sections on (0,7r) and 

(—71,0) and extend them to V^. We relabel them with Roman labels 
I and II as above. This is illustrated in Figure 2. Then equation (35) 
holds and in this basis the matrix of 0 becomes 

(36) 

«>„(k) = («(-k),e<(k))).,,=7,„ = -'=^"'■"'0 e f/(2) 

This matrix is skew-symmetric at fixed points k = —k and its Pfaffian 
there, which is just the upper right corner, is Pf(tc„(k)) = — 

This type of argument also works for other regions V = V~^UV~ with 
T{V^) = V~, Vo := V~^ n F“ C Fix{T) s.t. all the path components 
of F \ Fo are contractible and contained in either F"*" or F“, i.e. there 
is no path from F^ to V~ that does not go through a fixed point. 


3.2.3. Definition. In this subsection, we will briefly recall the definition 
of the Z2 invariant u using time reversal polarization following [32]. 
Originally a Z2 invariant was introduced by Kane and Mele [44] in the 
study of the quantum spin Hall effect in graphene, whose geometry is 
that of the honeycomb model with spin-orbit interaction. A simplified 
account was given in Fu-Kane [32]. The Kane-Mele invariant will be 
interpreted as a mod 2 analytical index, so that it counts the parity of 
the spectral flow of edge states in Z2 topological insulators [41]. 

Given the decomposition of H = 0 and restricting each Vn over 
an S'^, the above sections can be used to define an electric polar¬ 
ization in a cyclic adiabatic evolution. To this end, define the partial 
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polarizations as the integral of the Berry connections, 


/ TT 

al^{k)dk, s = IJI 


where the Berry connections are locally dehned as a* (/c) = i{u^{k)\dk\u^{k)). 
Using the combinations of the partial polarizations, one fnrther dehnes 
the charge polarization Pp and the time reversal polarization P 0 , 

Pp ■= pi + pii^ Pp ■= pi - pii 


For time reversal invariant topological insnlators, the charge polar¬ 
ization is always zero, and the time reversal polarization delivers the 
topological Z2 invariant. 

In particnlar, one compntes the time reversal polarization for the 
codimension 1 snbmanifold {(A;,0)|fc G [—tt, vr]} C in 2d. This 
corresponds to an edge (or edge states) of the effective Brillonin zone, 
which is a fnndamental domain for the time reversal Z2 action. This 
embedded has two hxed points, say F = {0, vr} C [—vr, tt]. One 
compntes, see e.g. chapter 4 of [63] : 


Pe = 




1 L det(wn(7r)) _ 2 Pf(w„(7r)) 

2Tri det(w„(0)) Pf(tc„(0)) 


Here Wn is the transition matrix (36), After exponentiating both sides 
and nsing the identity = —1, one obtains 


= TT V det(-u;n(rt)) 

^ “ 11 Pf(n;„(P.)) 

i —U, 7 r 


JJ s5(n(Pf(w„(Pi))) 

ri=0,7r 


It is important to point ont that det(wn{k)) is a continnons choice 
of sqnare root dehned on the interval from 0 to tt. This is possible 
by extending a choice of a branch of the sqnare root at 0 along the 
interval [0,7r] as in [32]. This is what gives the coherent choice of 
normalization mentioned above, which manifests itself as a choice of 
the sign of the Pfafhan. Since both the Pfafhan and the chosen branch 
are sqnare roots of the determinant, they differ only by a sign ±1. 
There is a more highbrow explanation in terms of sections which will 
be presented snbseqnently. 

In terms of the time reversal polarization, the Kane-Mele invariant 
on relative to a choice of sqnare root of the determinant is dehned 
by 


v = Pe (mod 2) 
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In general, the Z2-valued Kane-Mele invariant v is defined by 


(37) 


. Njy _ TT a/ det w„(rj) 

^ “ 11 Pf(in.(r,)) 


for the fixed points F of the time reversal symmetry. Here, again, one 
has to be careful to say how the branches of the square root of det are 
chosen. This is dictated by the embedding of S^s, which pass through 
the fixed points. In the 2d case these are chosen to be {(fc,0)} and 
{(fc,7r)} in [32]. This is actually not quite enough as realized in [32] 
and explained in [33]. Namely one needs a branch of the square root 
defined on all of T^. Since Ci{detB) = 0 there is indeed a global section. 
One can then argue that this section has a square root [31, 26], and 
below 3.3.3. Notice that this is possible, but it is not possible to find a 
basis of global sections, since the bundle need not be trivial. Unlike the 
case for for the 2d invariant, the result is actually invariant under 
such a choice of section, see [33, 31]. The 3d case is discussed in the 
next section. 


3.2.4. Spectral flow. In the appendix of [32], Fu and Kane showed that 
the 2d Kane-Mele invariant is equivalently given by 



where a and / are the Berry connection and curvature, EBZ is the 2d 
effective Brillouin torus with open cylindrical geometry and its bound¬ 
ary dEBZ consists of two circles. The equivalence can be proved by 
expressing the partial polarizations by integrals of the Berry connection 
and curvature. One has to calculate modulo two since this definition 
of V is not gauge invariant, indeed, the integral of a over the two circles 
always changes by an even number under a gauge transformation. 

An edge state explanation of the invariant is as follows: the 
can be viewed as chiral edge states. At any fixed point these cross 
each other and, after fixing the Fermi energy to be Ep = 0, change 
signs. By the time reversal symmetry, it would be enough to consider 
one chiral edge state, since the other one can be easily recovered by 
reflection. The Pfaffian would produce a negative sign if the spectral 
flow of the chosen chiral edge state is changed by -|-1 or — 1 under the 
gauge transformation Wn- Thus v keeps track of the parity of the change 
of signs of the Pfaffian under the time reversal symmetry represented 
by Wn, which can be interpreted as the evenness or oddness of the 
spectral flow or a mod 2 analytical index. This is made precise below. 
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In physical terms, a negative sign of the Pfaffian also corresponds to 
the existence of an unpaired Majorana zero mode. Under the bosoniza- 
tion procedure, two Majorana zero modes pair together to create an 
effective composite boson, so only the parity of Majorana zero modes 
matters. For z/ = 1, the effective theory consists of a collection of com¬ 
posite bosons (or none) and an unpaired Majorana zero mode such as 
a single Dirac cone. 

3.2.5. Strong Z2 invariant. We have dehned the Kane-Mele invariant 
or the analytical Z2 index above precisely in 2d. One can generalize it to 
3d [33] and higher dimensions using the formula (37), which is called the 
strong Z2 invariant. But, one can also consider various subtori and their 
Z2 invariants, which are called weak Z2 invariants. This in turn allows 
one to dehne the strong invariant rigorously via recursion. We will go 
through this in detail for a 3d topological insulator, whose momentum 
space (as the bulk space) is assumed to be with T(k) = — k. There 
are eight time reversal invariant points in total: (ni, 77.2,713) : Ui G 
{0,7r}. 

The three standard emdeddings of —)■ which have constant Tth 
coordinate equal 0,7r each contain 4 hxed points and give rise to six 
Z2 weak invariants ^'^^=0, r'ni=n,i = 1, 2, 3, computed by the restriction 
to the particular as a 2d invariant. These 2d Z2 invariants are not 
independent and the relation between them is given by 

^_X)^"2=0(_= ^_X)^"'3^°(_l)^"3“'''' =; (_1)^° 

which collects complete information from eight hxed points for different 
decompositions of into the strong invariant z/q. 

In practice, physicists use four Z2 invariants (z/q; P 1P2P3) to charac¬ 
terize a 3d topological insulator, where uq is the 3d Z2 invariant and 
z/j {i = 1, 2,3) are 2d Z2 invariants. 

For 3d physicists often use the momenta ki to he in [0, 2] mod 2 and 
the multiplicative form of Z2 instead of the additive form, i.e. u = ±1 
instead of z/ = 0,1. Following these notations (z/q; r'ir' 2 ^ 3 ) are, 

(39) Uq := Z/„.=o^'ni = l, ■= ^2 ■= J^n2=l ^3 ■= J^n3=l 

A product such as above indicates this multiplicative notation. 

For instance, a weak topological insulator could have z/i = 1, z/2 = 0, 
z/3 = 1 but z/q = 0, which is denoted by (0; 101). They are called weak, 
because they do not have to be stable under general deformations. We 
will talk about these four Z2 invariants again when we discuss the KR- 
groups of the 3-torus in section §5.4. 

The 3d Kane-Mele invariant also counts the parity of Dirac cones, 
i.e., conical singularities produced by surface states at the hxed points. 
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where we regard the codimension-1 sub-manifolds as edge spaces 
of the effective Brillouin zone. If there exist even number of Dirac 
cones, then the 3d Z 2 invariant is zero, i.e., vq = 0, but some 2d 
z/j {i = 1, 2, 3) could be non-zero when restricted to a submanifold T^. 
In this case, such time reversal invariant material is called a weak 
topological insulator. On the other hand, if there exist odd number of 
Dirac cones, then z/q = 1 and it is called a strong topological insulator. 

3.2.6. Holonomy. In this subsection, we will express the Kane-Mele 
invariant (originally in a quotient form) as a difference form with the 
help of holonomy. 

We will consider the Bloch bundle B restricted to R^, the effective 
boundary of the EBZ . In this region, there are no Kramers pairs, 
since 0 does not act hber-wise. But we can consider either |M^^(k)) or 
0 |m^(— k)) as interpolating between the quaternionic bases at 0 and vr, 
see §3.3.1. 

Dehne h as 

(40) h=-[ a^\k) - a\-k) dk = -(x„(7r) - Xn(0)) 

Jo tt 

with the help of the relation al^{k) — al^{—k) = dkXn{k) [32]. 

This can be either seen as a half quaternionic phase winding or a 
quaternionic holonomy around the loop —n —>■ vr. This holonomy 
basically gives the Pfaffian part of the time reversal polarization Pg 
discussed in the previous subsection. 

Let us start with h as the Pfaffian of the time reversal operator 0 and 
express the Kane-Mele invariant u as the difference of two geometric 
objects. Since Wn{7i) and Wn(0) are skew-symmetric matrices, their 
Pfaffians are 

Pf(u;„(7r)) = Pf(w„(0)) = 

Thus 

lnPf(w„(7r)) - lnPf(M;„(0)) = 2{k - i)Tri - f(Xn(7r) - Xn(0)) 

for some integers k,i & Z. Roughly, the holonomy is the difference of 
the logarithmic Pfaffians up to some winding numbers, 

vrh = i[ln Pf(r(;„(7r)) — In Pf(t(;„(0))] -|- 27r{k — i) 

To eliminate the ambiguity in choosing the integers from the complex 
argument of the logarithmic function, we modulo both sides by 27r, 

Pf(tnn(7r)) 

Pf(lTn(0)) 


7rh = i In 


(mod 27r) 
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It is easy to get an equality by exponentiating as before, 

PfK(o)) 


(41) 


- 1 )" = 


Pf(t(;„(7r)) 

The Pfaffian of Wn is well-defined only at the hxed points k = 0,7i, 
but the determinant of Wn is dehned everywhere for k G [0,7r]. Let us 
start with the equation satished by the determinant, 


1 

2 


(i[ln det(t(;„(A;))] = 2 


det(wn(7r)) 

det(wn(0)) 


= In 


A/det(wn(7r)) 

A/det(wn( 0 )) 


For any skew-symmetric matrix A, the sign of the Pfaffian satishes, 
FfA = det A, sgn(P{A) = Vdet A/PfA 
Then the above can be written in the Pfaffian and its sign. 


A/det(w„(7r)) ^ 


Pf(w„(7r)) 


In 


sgn(Pi{wn{7T))) 


v/det(wn(0)) LPf(«^n(0))J [s^n(Pf(M;„(0)))J 

Change the determinant into a trace, we have 


1 

2 


(i[lndet(tCn(/^))] = X / d[tr In Wn{k)] = - / dktr{w^ dkWn) 


Therefore, we have 
Pf(Wn(T)) 


In 


LPfK(o)) 


-|- In n sgn{Pi{wnia))) = - dktr{wJdkWn) 


a=0,7T 


In order to compare it with the holonomy, we exponentiate it. 


Pf(wn(7r)) 

Pf(lTn(0)) 


= (-l) 




JJ s5(n(Pf(M;„(ri))) 


ri=0,7r 

where the half winding number n is dehned by 

1 


(42) 


n(«'n) At / dktr{Wn'-dkWn) 

./n 


Then the above is equivalent to, 

(-l)-h = (-1)'^+^ 


In other words, the Z 2 invariant v is equivalently dehned by 
(43) z/ = n — h (mod 2) 

In contrast to the Kane-Mele invariant being a quotient in the previous 
subsection, this diherence form is another way to compare the square 
root of the determinant (as a half winding number) with the Pfaffian 
(as a holonomy). 
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3.3. Topological band theory. With the additional structure given 
by the time reversal symmetry, the Hilbert bundle modelling the band 
structure becomes a Quaternionic vector bundle. Specihcally, the time 
reversal Z 2 symmetry introduces an involution on the Brillouin torus 
and a Quaternionic structure on the physical Hilbert space, that is the 
sections of the Bloch bundle. Over the time reversal invariant points, 
the quaternionic structure gives rise to a symplectic structure of the 
restricted Bloch bundle. 

3.3.1. Quaternionic structure. In the bundle theoretic framework, one 
obtains the Hilbert space as the space of sections of the Bloch bundle 
Ti = r(i3). Since the bundle B splits as 0„Vn, so does the Hilbert 
space as with Hn = r(Vn). 

With the extra structure of the time reversal operator 0, the complex 
Hilbert space 'H(C) = "H can be viewed as a quaternionic Hilbert space 
^(EI), by setting j = 0, with = 0^ = —1. The action of El = 
span^QQ, j, k) = spanc{l,j) is given by i, i.e. multiplication by i, j 
acting as 0, which is anti-unitary and hence anti-commutes with i, 
and k = ij = iQ. 0 respects the decomposition of "H, so that each Tin 
becomes a El vector space 'Hn(EI) and 'H(EI) = 0„'Hn(EI). The hbers 
of Vn over a hxed k are 2-dimensional complex vector sub-spaces. 0 
in general permutes these sub-spaces, except over the hxed points. At 
these points the hber becomes a quaternionic vector space. That is 
B\t is a quaternionic bundle and each hber VnVi — a 2d C 

vector space becomes a Id El vector space. In particular "Hr = r(H, F) 
is an El invariant subspace. Over F, if \un) is a (local) basis over El 
then \un) and \Qun) is a basis over C. These are usually put into a 
so-called Kramers pair 

l^n) := (|Mn),0|Wn)) G Hr(H) C HniM) 

then picking {\un), Olwn)) as a basis, we see that 
0 ■ = JT„, J = i(Jy = 

where is viewed as a 2-vector. 

The relation to the Pfafhan is as follows: dehne a bilinear form on 
the complex Hilbert space "H by 

(44) a;0(0,i/>) = (00|i/>) 

This form becomes a symplectic form for Vn over the time reversal 
invariant points; ca© G A^(Vn)*. Thus Vn|r —)■ F is a symplectic vector 
bundle. The associated almost complex structure is Jn = —iay. Hence 
the Bloch bundle restricted to the hxed points Vn|r F carries an 


0 1 

-1 0 
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associated almost complex structure Jq = 0^ Jn- We have that 0 = 
KJ. 

Recall that the compact symplectic group Spin) is defined as the 
quaternionic unitary group, 

M e Sp{n) = U{n, H), he., MM^ = M = 1 

where = irflji) is the combination of taking the matrix trans¬ 

pose and the quaternionic conjugate. In particular, one has S'p(l) = 
SU{2). The Berry connection, which has the local form an(k) = 
(T„(k)|d|T„(k)), satisfies the identity = 0, i.e., G sp(l). 

Remark 1. Instead of (44) matrix coefficients of Q as a morphism 
T-i ^ if are 

Cu((/),'0) = ((/>,0i/’) = (0|0|V') 

This is actually an anti-linear 2-form and a section of It is 

related to ojq by 

(45) 

3.3.2. Quaternionic K-theory. The extra structure that 0 adds to the 
bundle has actually been studied by mathematicians and is called a 
Quaternionic structure [28]. In this setup one considers a locally com¬ 
pact topological space X together with an involution t. The relevant 
example here is the time reversal symmetry acting on the Brillouin 
torus. A Quaternionic bundle on such a space is a bundle B together 
with an anti-linear morphism y, which is compatible with t and satis¬ 
fies = — 1. In our case y = 0 and the Quaternionic bundle is simply 
(S,0). 

If instead y^ = 1, then the pair is called a Real bundle, which was 
first considered by Atiyah [4]. The capital “Q” and “R” indicate that 
the morphism y does not necessarily act as a bundle isomorphism, but 
rather interchanges the fibers over x and t{x). If the involution y is 
trivial, one obtains quaternionic and real bundles. Looking at classes 
of bundles one obtains the respective iL-theories. Complex K = KU, 
real KO, quaternionic or symplectic KSp, Real KR and Quaternionic 
KH (or KQ). These are alternatively defined using homotopy theory, 
the Grothendieck construction or stable equivalence [4], see below. As 
we consider the case with 0^ = —1, we will be mainly interested in 
Quaternionic structures. The others appear when there are other sym¬ 
metries present, see §5.1 below. For instance the case with 0^ = 1 
represents particles with integer spin [71]. 

Denote KH{X) the K-group of Quaternionic vector bundles over 
an involutive space (X, r). Sometimes the notation KQ{X,t) is also 
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used [26]. If the involution r is trivial, then KH(X) is reduced to the 
quaternionic K-theory KSp{X) = Km{X) = KQ{X,idx)- Similar to 
KR-theory, the KH-groups have a period 8, 

KH\X) ~ KW^^{X) 

The relation between the KH-groups and the KR-groups is given by 

KH\X) ~ KR-^{X) 

Since the Bloch bundle (H, 0) —)■ (T'^, T) has a Quaternionic structure 
(0^ = —1) rather than a Real structure, it is better to use KH{T’^) 
for time reversal invariant topological insulators. However, the above 
relation tells us that if we shift the KR-groups by 4, then the KR- 
classihcation [31] agrees with the KH-classi£cation [26] for the clas- 
sihcation problem of Bloch bundles over the Brillouin torus for time 
reversal invariant models. 


3.3.3. Determinant line bundle. In this subsection, we hrst give an¬ 
other equivalent dehnition of the Kane-Mele invariant by looking at 
the determinant line bundle following [31]. Then we discuss the closely 
related FKMM-invariant in equivariant cohomology following [26] . The 
calculations of §3.2.6 yield the equivalence to the classical dehnition of 
the Z 2 invariant using integrals. 

For the rank 2N Bloch bundle {B, 0) —?• T'^ consider its determinant 
line bundle, which is given by: DetB := R^B —>■ T'^. Taking the 2N- 
th exterior power DetQ of 0, we obtain an anti-linear morphism that 
squares to -|-1. This yields a Real bundle {Det B, DetQ), and when 
restricted to the hxed point set F a real structure {Det B\t, DetQ\r). 
The DetQ invariants of DetB\r form a real sub-bundle, which is de¬ 
noted by Det^Br- 

This real bundle has an orientation given by the hrst component of 
the Bloch states. This means that if the \un) are a quaternionic basis, 
then |mi) A 0|mi) A ■ ■ ■ A \un) A 0|MAf) is a canonical orientation. It is 
Det 0 invariant, since Det 0(0 A 00) = 00 A0^0 = —00A0 = 0A00. 
It is canonical since the space of quaternionic bases is connected. 

Alternatively, the symplectic structure uq defines an anti-linear 2N 
form —aka. the Pfafhan form— which is 


Pf(a;0) := 


, ,AN 

N\ 


uq a ■ ■ ■ a ca© 

iv! 


It restricts to an M-linear form on the real bundle DeteBr yielding an 
orientation. This follows from ue{Qu, Qv) = we{u,v). 

As mentioned above, Ci{DetB) = 0 and hence there is a nowhere 
vanishing section. We can make this section real, i.e. land in the Det 0 
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invariant part by averaging, see e.g. [31] for details. We will denote 
this section by s. By comparing s(rj) with the canonical orientation 
of the real line Det^Br, one dehnes the Z 2 invariant as the product of 
signs of the invariant section 


u = 


n =n 

Uer uer 


Pf(ri) 


where Pf : F —)■ Det^Br is a global section representing the canoni¬ 
cal orientation determined by the Pfaffian form. he. the sign is posi¬ 
tive/negative if s is in the positive/negative half line determined by the 
orientation. Abusing notation, one can set s = \/det since it is DetQ 
invariant, and rediscover the Kane-Mele invariant. In other words, the 
Z 2 invariant can be dehned equivalently by comparing the orientations 
of the restricted determinant line bundle and the Pfaffian line bundle 
over the fixed points. 

In [26], a so-called FKMM-invariant was constructed almost along 
the same lines as above. The restricted determinant line bundle DetBlr 
is trivial since F is a finite set. Furthermore there is a canonical iso¬ 
morphism 

detp : Dct B\y — F x C 

This section is related to the above equivariant section by \/det = 
The Real determinant line bundle (Det B, DetQ) has an associ¬ 
ated equivariant trivialization hdet '■ DetB = x C, and its restriction 
hdetir gives another trivialization of DetB\T- The composition 

hdet\T ° : F X C —)■ F x C 


defines a map 

(46) cjb : F-)■ Z 2 , s.t. {hdet\T ° detY^){x,X) = {x,ub{x)\) 


Then the FKMM-invariant is defined as the equivariant homotopy class 
of ub up to equivariant gauge transformations 

(47) k{B) := K) e [r, U(l)\zJ[T\ t/)!)]*. 


In general, k defines an injective map from Quaternionic vector bun¬ 
dles to relative equivariant Borel cohomology classes over an involutive 
space (X, <^), 

(48) « : Kecf2-(X,<^) ^ Z(l)) 

See §4 of [26] for more details about the injectivity, the equivalence 
between the FKMM-invariant and the Kane-Mele invariant was also 
proved in that section. 
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3.4. Homotopy theory. The study of homotopy theory of time re¬ 
versal invariant Hamiltonians can be traced back to the 80’s [12, 13]. 
Denote the space of non-degenerate Hamiltonians, i.e. non-degenerate 
k X k Hermitian matrices by Sjk- It is homotopic to U{k)/U{1)^. The 
lower homotopy groups were computed in [13], 

vri(i3fc) = 0, TT2{S)k) = 

In particular, for the 2d case (i.e., k = 2) 712 (^ 02 ) = which reflects 
the existence of the TKNN integers (or the hrst Chern number) as the 
Z-invariant. 

As mentioned previously, from [68], we know that for a family of 
Hamiltonians H G io, its eigenvalue degeneracy is generically codimen¬ 
sion three, so in 3d a 8“^ can be used to enclose the monopole and the 
hrst Chern number of the restriction of the Bloch bundle provides a 
topological invariant, see §2.3. This corresponds to the fact 712 (i 02 ) = Z- 

If one includes time reversal symmetry, two more spaces become of 
interest. The hrst is which is the subspace of Hermitian matri¬ 
ces allowing possible pair degeneracies, that is Eigenvalues 1 and 2 or 
Eigenvalues 3 and 4 may be degenerate and so forth. This is the space 
for families of Hamiltonians which have at most Kramers degeneracy. 
The second is the subspace 

(49) = {H-1 QHQ-^ = H} 

where now 0 is the given action of time reversal on the Hilbert space, 
that is an anti-unitary operator with 0^ = — 1. To not overburden the 
notation we omitted the index A;, but we will always assume that the 
dimension of the family of Hamiltonians is hxed. Moreover, for time 
reversal invariant systems the interesting case is when the Bloch bundle 
has even rank, so we will assume that k = 2n. In this case, independent 
of the particular choice 0, each element in i^® can be represented by 
a non-degenerate quaternionic n x n Hermitian matrix, and the space 
i^® is homotopic to Sp{n)/Sp{l)'^, see e.g. [12]. 

The relation is that the family of Hamiltonians over the Brillouin 
torus can be viewed as a map if : T'^ —)■ If we add the time 

reversal symmetry, we additionally have the two Z 2 actions T on 
and conjugation by 0 on i^. The map H then becomes a Z 2 equivariant 
map. The condition that there are only the degeneracies forced by 
Kramers degeneracy implies that in fact H : S^Kr- Since this 

map is still equivariant it is determined by the restriction H\ebz '■ 
EBZ —)■ S^Kr- Said differently, under the time reversal symmetry, we 
only consider the half Brillouin torus giving a fundamental domain for 
the Z 2 action, i.e. the EBZ, since the other half can be easily recovered 
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by reflection. Furthermore, at the hxed points Fj G F the Hamiltonians 
are invariant, so that H\y ^ C Sj. 

The lower homotopy groups of time reversal invariant Dirac Hamil¬ 
tonians are known, see [12]: 

= 7r2(-^®) = = 0, M^n) - 

Analogously to the above, for a generic time reversal invariant family 
of Dirac Hamiltonians H G iof, its eigenvalue degeneracy has codimen¬ 
sion hve and Kramers doublets are distinguishable in 4d families, e.g. by 
enclosing the degeneracies by 5^. In other words, the homotopy classes 
of maps [S'^ , = 7r4(i3f) = Z are labeled by Z, which is just the 

second Chern number C 2 of the Bloch bundle. Thus C 2 manifests itself 
as a fundamental invariant for time reversal invariant models, which 
follows from 7r4(ij®) = Z. As before this yields no direct invariants in 
dimensions less than 4. For other roles of S''^, see §3.1.2. 

In [55] , Moore and Balents proposed that one could cap the two ends 
of the cylinder C = x I, representing the effective Brillouin zone to 
obtain a topological invariant. This capping turns C into an S'^ and 
one can then regard an induced family of Hamiltonians on 8“^, which 
however depends on the chosen extension to the caps. By a homotopy 
argument which will be explained below, they obtained the topological 
Z 2 invariant. 

First, Moore and Balents assert that vii= 0, which implies that 
the map H restricted to any of the two circles in the boundary of the 
effective Brillouin zone to is homotopic to the constant map. It is 
claimed in [55] that this homotopy can be made into a homotopy with 
two additional properties, which is what they dehne to be a contraction. 
The hrst is that the homotopy is actually a homotopy of equivariant 
maps, and secondly, that homotopy f{9, A) : x J —)■ S^xr satishes 

that f{6,'n') G This can be made rigorous by considering relative 
homotopy groups ,S}Kr, Hq) and the corresponding long exact 

sequence, see e.g. §4.1 in [42]: Consider H\si —)■ S^xr where is 
one of the boundaries of the EBZ cylinder. This map is actually still 
Z 2 equivariant, so that it is hxed by its restriction to the upper half 
circle, which is an interval. The restriction to this interval is a path 
p from Hq G to a point H.„ G and hence in 7ii{S}^,Sjxr,Ho). 
From the long exact sequence of relative homotopy groups this group 
is trivial and hence the path can be homotoped through paths pt to 
a path Pi entirely inside Considering the loops k = pt for points 
in the upper half-circle and It = QptQ~^ for the lower half-circle, we 
obtain an equivariant deformation of H\si to a map /i : —)■ with 

/t(7r) G Since 7ri(i5®) = 0, we can further contract to the constant 
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loop, again in an equivariant fashion, obtaining, finally, a contraction 
in the above sense. 

Hence the map from the cylinder representing the effective Brillouin 
zone can be replaced by a map from the 2-sphere, obtained by con¬ 
tracting the boundary circles to a point. More precisely, any map 
H : EBZ —)■ S^Kr is homotopic to a map that factors through the 
quotient space = EBZ/ ~ where ~ is the equivalence relation 
identifying all points in each boundary circle. This can also be visu¬ 
alized by gluing two caps onto the cylinder. The resulting map from 
H : ^ S^Kr produces a hrst Chern number, which will be denoted 

by Ci{H). This hrst Chern number depends on the choice of homotopy, 
though. What is argued in [55], is that choosing different homotopies 
changes the Chern number by an even integer: Ci{EI) = Ci{EI') mod 2. 
The argument relies on the fact that homotopies can be chosen to be 
contractions and that the lower homotopy groups tti, 712 of vanish. 

After eliminating the ambiguity of extra even Chern numbers from 
choosing contractions, a topological Z 2 invariant is obtained by capping 
the cylinder. This procedure to eliminate the ambiguity by a homotopy 
argument is in line with the solution to the gauge ambiguity in (38). 

4. Z 2 INVARIANT AS AN INDEX THEOREM 

In this section, we show how the Z 2 invariant can be viewed both 
as an analytical index based on the Maslov index and as a topological 
index by the Chern-Simons invariant. We have explained that the 
Kane-Mele invariant counts the parity of the spectral flow of edge 
states in §3.2.4 and that this in turn, in the language of index theory, 
is the analytical Z 2 index. The topological Z 2 index is the focus of 
the second subsection. In physics, it corresponds to the Chern-Simons 
invariant, more precisely, the mod 2 version of the topological term 
in the Wess-Zumino-Witten (WZW) model. The last subsection is a 
summary of the relations between the different interpretations of the 
topological Z 2 invariant. 

4.1. Analytical Z 2 index. The mod 2 analytical index was hrst intro¬ 
duced by Atiyah and Singer [8, 9] for real skew-adjoint elliptic operators 

P. 

(50) ind 2 {P) = dimker(P) (mod 2) 

which can be equivalently described by a mod 2 spectral how. We 
argue that the Kane-Mele invariant can be interpreted as a mod 2 
analytical index of the ehective Hamiltonian of Majorana zero modes, 
which is explained clearly in our work [47] . For the original Kane-Mele 
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invariant, a possible connection was first observed by Tony Pantev, see 
reference [14] in [44], 

In the quantum spin Hall system, the Kane-Mele invariant i/ = 0 is 
equivalent to the fact that the Fermi level intersects a chiral edge state 
even times, while z/ = 1 means the intersection number is odd. This is 
exactly the dehnition of the mod 2 analytical index (50). Here we put 
it into the symplectic setting and give an equivalent dehnition of the 
Kane-Mele invariant using the Maslov index as in [10], which gives a 
geometric realization of the analytical index. 

4.1.1. Maslov index. Assume that the Fermi level sits inside the band 
gap between the conduction and valence bands, in this subsection we 
view the Fermi level as a hxed reference Lagrangian subspace M C 
From the helical edge state, we pick one chiral edge state and model 
it by a continuous piecewise-smooth real-valued function, which is also 
modeled by one Lagrangian subspace or several Lagrangians in If 
we count the intersection points of the chosen chiral edge state with 
the Fermi level, then the 2d Z 2 invariant can be equivalently dehned 
by the Lagrangian intersection number or Maslov index as an edge Z 2 
index [10]. 

It is well known that the Maslov index can be characterized by a 
set of axioms, which has different geometric and analytical realizations 
such as eta invariant or spectral flow. In this subsection, we recall 
the relations between the Maslov index, eta invariant and spectral flow 
within following [21]. More generally, on manifolds with boundary 
the spectral flow of a family of Dirac operators can also be computed by 
the Maslov triple index, the interested reader may consult for example 
150] . 

In the simplest symplectic vector space with the standard 

symplectic form cu, L C is a Lagrangian subspace, or simply La¬ 
grangian, of if L = L-*“, which is equivalent to L being 1- 

dimensional in Denote the 1-dimensional Lagrangian Grassman- 
nian Gr{l, 2) in (M^, uj) by A(l), and each element in A(l) is represented 
by a line L{6) with slope 9. The square function 

Sq : A(l) ^ 5^ L{9) ^ 

induces a diffeomorphism A(l) ~ U{1)/0{1) — and the funda¬ 
mental group of A(l) is 7ri(A(l)) = = Z. 

In (M^,a;), the geometric Maslov index is dehned as the intersection 
number of a line L{9) G A(l) with a reference line usually hxed as 
M = L(0). Remember that in the 2d quantum spin Hall system, we 
model the chiral edge states by Id Lagrangians, and we are interested 
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in their intersections with the Fermi level. One version of analytical 
Maslov index is defined by the eta invariant of a real self-adjoint Dirac 
operator D{Li,L 2 ) where Li,L 2 G A(l). Another version of analyti¬ 
cal Maslov index is defined by the spectral flow of a family of Dirac 
operators {D{Li(t), L 2 (t)),a < t < b}. In other words, the geomet¬ 
ric Maslov index can be interpreted as the spectral flow of a family of 
Dirac operators, which is intimately related to the eta invariant. 

We now show how to compute the analytical Maslov indices explicitly 
in (M^,a;). Denote the unit interval by / = [0,1], let Li, L 2 ) be 

the Sobolev space as the completion of the space of paths connecting 
two Lagrangians Li,L 2 , i.e. 

0 e C\I, M^) S.t. 0(0) e Ai, 0(1) e L 2 

The Dirac operator D{Li,L 2 ) dehned by 

(51) D(Ai, L 2 ) : Li, L 2 ) ^ L\I, R^)- 0 ^ 

is a Fredholm operator and its kernel is given by constant functions 

Ker{D{Li, L 2 )) = {0(t) = constant point G Ai fl L 2 } 

When Re{s) 3> 0, the ? 7 -function of a self-adjoint elliptic operator D 
is dehned by 

( 52 ) = 

A>0 A<0 ^ ' 

which admits a meromorphic extension to the complex plane and is 
regular at s = 0. One dehnes the eta invariant as 77 (D) := r]{D, 0), and 
the reduced eta invariant 77 (D) as 

(53) f,(D) := ’KJ)+dim(kerJ) 

Roughly, the eta invariant is the number of positive eigenvalues minus 
the number of negative eigenvalues. In symplectic geometry, one dehnes 
the analytical Maslov index by the eta invariant 

(54) /Xa := 77 (D(Li, A 2 )) 

We stress that this analytical Maslov index is real-valued. 

For example, consider two typical Lagrangians M and L{6) = 
in (M^,a;), the path connecting them can be taken to be 

0 (t) = for some ieZ 



NOTES ON TOPOLOGICAL INSULATORS 


37 


where the integer i is added since the Lagrangian is invariant 

by multiplying a minus sign. Since 

D(R, L(e))(j) = = Utt + eu 

at 

the eta function of D{R, L{6)) is 

ri{D{R, L{9)), s) = ^ [Ch{s, 6 / 71 ) - Ch{s, 1 - O/tt)] 

where (h{s,x) is the Hurwitz zeta function (h{s,x) = 
for 0 < T < 1. The Hurwitz zeta function has a regular point at s = 0, 
indeed Cj7(0,a;) = ^ — x. Hence the real-valued Maslov index is easily 
obtained by direct computation 

(55) iaa = v{D{R,L{e)),0) = l-- 

71 

In order to be compatible with the integral Lagrangian intersection 
number, i.e., the geometric Maslov index, the integer-valued analytical 
Maslov index is defined as follows. Dehne the canonical 1-form on A(l) 
by pulling back the standard 1-form on S^, 


n := 



which is the generator of the cohomology group i/^(A(l), Z) ~ Z) ~ 

Z. Now the integer-valued analytical Maslov index is dehned for a pair 
of continuous and piecewise smooth Lagrangians f(t) = {Li(t), L 2 (t)), 
for some interval a <t < b, 


(56) /(„(/) fma - L;f!] + ^fAD(L„L,yt 

where 



-'Y][dsi{b) - 9si{a)], s = l,2 

71 —' 


In (]R^,a;), let Li = M, L 2 = then it is easy to compute the 

eta invariant 



1 

1 

for - f < 0 < 0 

v{D{L,,L2)) = 

{ 1-f 

for 0 < 0 < f 


1 0 

for 6^ = 0 

the pull-backs are given 
(M, Me**), it follows that 

by = 0, 

L 2 H = 99/77. For f{t) 

Aa(/l[-^,0]) = 0, 

ba(/|[0, J) 

= M/I[-^,5) = i 
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The spectral flow of a family of self-adjoint elliptic operators, in 
practice Dirac operators, Dt = D{Li(t), L2{t)),a < t < b is defined as 
usual by the net number of eigenvalues changing sign (-1-1 for increasing, 
— 1 for decreasing) from a to b, denoted by sf{Dt). A second integer¬ 
valued analytical Maslov index is defined by the spectral flow of the 
family of Dirac operators associated with the pair f(t) = {Li(t), L2(t)), 

(57) fiaif) ■= sf{D{Li{t), L2{t))), a<t <b 

Recall that D(]R, Me*®) has eigenvalues 7 ri + 6 for —j < 6 < When 
6 approaches 0 from the left, the eta invariant has the limit 


Meanwhile, we also have the right limit 


Immediately, the spectral flow is easily obtained from the jump of the 
eta invariant crossing 0, i.e. 

s/(D(Li(f), L2m = ^[r/(D(Li, L 2 ), 0+) - r]{D{L,, L 2 ), 0-)] = 1 

The analytical Maslov index unifies the spectral flow and the eta in¬ 
variant in the symplectic setting, that is, the spectral flow of a family 
of Dirac operators can be computed as the variation of eta invariants. 

To end this subsection, we define the edge Z 2 index [10] by the Maslov 
index /i modulo two, 

(58) 1/ ^ 5^ /i(M, L(6*j)) (mod 2) 

i 

where the Fermi level is fixed as M = T(0) and one chiral edge state 
is assumed to be piecewisely modeled by a finite collection of Id La- 
grangian subspaces L{di) G A(l). The equivalence with the Kane-Mele 
invariant can be seen by expressing the edge Z 2 index as the winding 
number of a closed unitary path, for more details see [10]. 

4.2. Topological Z 2 index. In this subsection, we will discuss how to 
compute the Z 2 invariant of a 3d topological insulator by the topolog¬ 
ical Z 2 index. One starts out with a 4d time reversal invariant model 
and applies dimensional reduction to obtain a 3d Chern-Simons held 
theory. From the 3d bulk Chern-Simons (CS) theory, one obtains a 
2d topological Wess-Zumino-Witten (WZW) term by the CS/WZW 
duality. It is shown in [31, 47, 69] that for a 3d topological insulator, 
the Chern-Simons invariant, i.e., the WZW term, and the Kane-Mele 


/ 26 

lim riiDiLi, L 2 )) = lim (1- 

6»-s>0+ e^O- \ 71 


( 26 

lim riiDiLi, L 2 )) = lim ( —1- 

0^.0- ^ ^ 0-5>O- y TT 
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invariant are equivalent, and for this reason the topological Z 2 invariant 
will be called the Chern-Simons invariant in this subsection. In string 
theory, the topological WZW term naturally appears in the bosoniza- 
tion process [73], and its global SU{2) anomaly explains why it is a 
mod 2 index [72] . 

Mathematically, the Chern-Simons invariant is computed by the odd 
topological index, that is, the integral of the odd Chern character of 
a specihc gauge transformation, which is an application of the odd 
index theorem for Toeplitz operators on compact manifolds. So the 
topological Z 2 invariant is identihed as the mod 2 reduction of the odd 
index theorem in 3d. Or equivalently, the Z 2 invariant is the mod 
2 version of the Witten index Tr(—1)^, where F counts the number 
of Majorana zero modes such as Dirac cones. In a modern language, 
the topological Z 2 index is obtained by pairing the K-homology of 
Majorana zero modes with the K-theory of the Bloch bundle, for details 
see [47]. 

4.2.1. CS/WZW correspondence. Zhang et ah [61] hrst considered a 
time reversal invariant (4+l)-dimensional model generalizing the (2+1) 
quantum Hall system, the effective action functional was obtained by 
linear response theory, 

(59) = 

where C 2 is the second Chern number 

(60) C 2 (f) = ^/ d^ke‘>“tr{fyft,) 

with non-abelian local Berry connection and curvature dehned as be¬ 
fore, 

tti = fij = diOj - djOi + i[ai, aj] 

In order to get the effective action functional of a 3d topological insu¬ 
lator, dimensional reduction was applied to the above effective action 
5*4^ and the second Chern number as well [61]. The idea of dimen¬ 
sional reduction is to replace one spatial or momental dimension by an 
adiabatic parameter, and then integrate out this adiabatic parameter 
carefully. 

For the external gauge held, one can reduce the dimension by hxing 
a convenient gauge such as the Landau gauge. As for the 2nd Chern 
number C 2 (/), one considers the Chern-Simons class because of its 
relation with the second Chern character on the level of characteristic 
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classes. Indeed, the Chern-Simons form 053 ( 0 , /) is the boundary term 
of the 2nd Chern character, for example see [56], 

d 053(0,/) = ch 2 {f) 

In addition, the Chern character can be expressed in Chern classes, in 
our case 

ch2if) = ^[cjif) - 2C2(/)] = -02(/) 

since oi(/) = 0 for time reversal invariant models. 

After dimensional reduction, the effective action for a 3d topological 
insulator is given by 

( 61 ) 

where the magneto-electric polarization P 3 is dehned as the Chern- 
Simons action of the Berry connection, 

Pd,{ai) = j (fkpPtr{aifjk - ^Oi[oj,Ofc]) 

It is related to the d-parameter known in condensed matter physics by 

(62) 9 := 27rP3 G { 0 , 7 r} mod 27r 

In general, the Chern-Simons action is gauge invariant up to a wind¬ 
ing number. More precisely, under a gauge transformation, 

a, I—)■ af = g~^aig — g~^dig, for g :T^ ^ U{n) 

one has 

(63) AP3 = Psi^af) -Psitti) = J d^ke^Ptr{g~^digg~^djgg~^dkg) 

the right hand side is an integer, namely the winding number of g. It 
is also known as the topological WZW term. 

The Chern-Simons invariant i; of a time reversal invariant system is 
dehned in [61] by the change of magneto-electric polarization under the 
gauge transformation induced by the time reversal symmetry modulo 
two. More precisely, if w is the specihc gauge transformation induced 
by the time reversal symmetry given by the transition matrix (36), the 
Chern-Simons invariant is dehned by 

(64) V = I [ d^ktr{w~^dw)^ (mod 2) 

2471"* Jj3 

Now the Bloch bundle splits as before, so that we have a block form 
of tc G P(2) X • • • X P(2) C U{2N). Furthermore, by the property of 
the transition matrix, the gauge group can be reduced from P( 2 ) to 
SU(2), see [69], hence v is the mod 2 version of the topological term 
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in the SU{2) WZW model. So the chiral edge states can be modeled 
by some current algebra, which is a Kac-Moody algebra in conformal 
field theory, the interested reader may consult [73]. 

Up to now, we considered the 3d Chern-Simons theory and the asso¬ 
ciated topological WZW term. A similar argument can be applied to 
the Id Chern-Simons theory, which can be viewed as a boundary theory 
of a 2d bulk theory. The time reversal symmetry also induces a specific 
gauge transformation, which can be plugged into the Id WZW term, 
i.e., the winding number. The mod 2 version of this Id WZW model 
can be used to study the continuous model of time reversal invariant 
Majorana chains. 

The physical picture behind the topological WZW term is the bosoniza- 
tion process, that is, two Majorana zero modes, such as two Dirac 
cones, pair together to become an effective composite boson. The topo¬ 
logical WZW term as an action functional describes a bosonic theory 
equivalent to the fermionic theory of Majorana zero modes. The built- 
in global SU(2) anomaly of the WZW term explains why we have to 
modulo two, i.e., the parity anomaly, which will be discussed in a later 
subsection. 

To sum it up, the CS/WZW duality gives the bulk-edge correspon¬ 
dence in 3d topological insulators. More precisely, the Chern-Simons 
action of the Berry connection describes a 3d topological insulator in 
the bulk, while the topological WZW term characterizing the edge 
states is the action of the specihc gauge transformation induced by the 
time reversal symmetry. In other words, the bulk field theory is the 
3d Chern-Simons field theory, and the boundary field theory is the 2d 
Wess-Zumino-Witten model. 

4.2.2. Odd Chern character. From the last subsection, we know that 
the Chern-Simons invariant can be dehned by the mod 2 WZW term of 
a specific gauge transformation. It is well-known that the WZW term 
gives rise to a topological index of Toeplitz operators. In this subsec¬ 
tion, we will review the odd Chern character of gauge transformations 
and its spectral flow following [37], so that the Chern-Simons invariant 
will be naturally interpreted as the mod 2 spectral flow through odd 
Chern character. 

For two connections Mo and Ai on some vector bundle, the relative 
Chern-Simons form is dehned by [37] 


( 65 ) 
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where 

A-t = (1 ~ + tvAi, A-t = Ai — *4-0 

There exists a transgression formnla connecting the Chern-Simons 
form and Chern characters, 

(66) dcs{AQ,Ai) = ch{Ai) — ch(4.o) 

where the Chern character is defined as nsnal, ch{A) = tr{e-^^). 

Fnrthermore, given a map from a compact {2k + l)-dimensional spin 
manifold M to the nnitary gronp g ■. M ^ U{n), its homotopy class 
is an element in the odd iF-gronp, i.e., [g] G K~^{M). The odd Chern 
character of g is dehned by [37] 

A 

(67) ch{g) := ^ ^y tr[{g-^ dgf^+^] 

which is a closed form of odd degree. Or eqnivalently, ch{g) is the 
relative Chern-Simons form cs{d, g~^dg) by Taylor expansion. As a 
classical example, the degree of g : —)■ U (n) is given by 


deg(^) 



ch{g) 


Replacing connections by Dirac operators, the role of the relative Chern- 
Simons form cs{d, g~^dg) is played by the spectral flow sf{D,g~^Dg) 

137], 

The analytic spectral flow of a Dirac operator D on M can be com- 
pnted as 

( 68 ) sf{D,g~^Dg) = ^f tr{Dte~^^)dt 

V ^ Jo 

where 

Dt = {I-t)D+ tg-^Dg, Dt = g~^[D,g] 

The Dirac operator D dehnes a Fredholm modnle in K-homology, and 
the spectral flow compntes the Fredholm index of the Toeplitz operator 
PgP by the pairing between K-homology and K-theory [37], 

(69) index{PgP) = {[D], [g]) = -sf{D,g-^Dg) 

where P := (1-1- D\D\~^)/2 is the spectral projection. Banm and 
Donglas [15] first noticed the odd Toeplitz index theorem, which is an 
identity connecting the analytical index and topological index. 


sf{D,g ^Dg)= f A{M) A ch{g) 

J M 


( 70 ) 
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where A is the A-roof genus, for a generalized odd index theorem for 
manifolds with boundary see [24]. 

In particular, we have A(T^) = 1 since A is a multiplicative genus 
and A(S^) = 1 for spheres. Hence the degree of g can be computed as 
the spectral flow on the 3d Brillouin torus, 

(71) sf{D,g-^Dg) = - ^ 

Finally, we apply the odd index theorem (71) to the Dirac Hamiltonian 
77 of a 3d topological insulator and the transition matrix w (36), then 
the Chern-Simons invariant is identihed with the spectral flow modulo 
two, 

(72) V = sf{H,w~^Hw) mod 2 

For manifolds with boundary, a well-defined Dirac operator involves 
the APS boundary condition, see [6], [50]. The analytic r^-invariant of 
a Dirac operator D is dehned by 

(73) ',(7)):= A 

A Jo vs 

If Dp{t) is a family of Dirac operators dehning a good boundary value 
problem, P is a spectral projection, then the spectral flow of Dp{t) can 
be computed as variations of the eta invariant, see [37], [50], 

sf{Dp{t)) = fj{Dp{t))\l 

where fj is the reduced eta invariant. As a remark, if a superconnection 
A is used in (73) instead of P, then it would define the eta form, which 
has important applications in the family index theorem and Witten’s 
holonomy theorem [19]. 

4.2.3. Green’s function. In this subsection, we define the Chern-Simons 
invariant of time reversal invariant topological insulators by Green’s 
functions following [70]. This approach has the advantage that such 
topological invariant is still valid even in interacting or disordered sys¬ 
tems. The idea of applying Green’s functions to the study of domain 
walls in D-branes or edge states in topological insulators can be traced 
back to early works by Volovik, for example see [67]. 
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By evaluating one loop Feynman diagrams, the second Chern class 
(60) in the 4d action (59) is equivalently dehned by fermionic propa¬ 
gators [38], i.e.. Green’s functions, 

(74) 

2 p 1 1 

C2 = — —^ef^^>’^Hr[GdM-^Gd,G-^GdpG-^Gd^G-^GdrG-^] 

15 J [271^ 

where the imaginary-time single-particle Green’s function of the Hamil¬ 
tonian H is dehned as 

G{iu,'k) := {iu — H{k))~^ 

For 3d time reversal invariant topological insulators, the electro¬ 
magnetic polarization P 3 is again obtained by dimensional reduction, 

(75) 

P 3 -=qJ du j -j^e^^'^P^tr[GdpG-^Gd^G-^GdpG-^Gd^G-^Gd^G-^] 

By a similar argument as the WZW extension problem, which will 
be explained in the next subsection, 2 P 3 is well dehned up to certain 
integer belonging to the homotopy group 7 r 5 (GL(n, C)) ~ Z. In order 
to get the topological Z 2 invariant in 2d, we need two WZW extension 
parameters and to reduce one more dimension. The 2d topological 
order parameter is dehned similarly, 

(76) 

P 2 ;= — / dudv / —^e^^''nr[GdpG-^Gd^G-^GdpG-^GduG-^Gd^G-^] 
120 J {271^ 

The same logic shows that P 2 induces the topological Z 2 invariant 
and P 2 = 1/2 corresponds to non-trivial topological insulators in 2 d. 

4.2.4. Parity anomaly. The parity anomaly caused by unpaired Majo- 
rana zero modes hts into an SU{2) gauge theory when the time reversal 
symmetry is translated into a specihc gauge transformation as discussed 
in §4.2.1. Fukui, Fujiwara and Hatsugai [34] explained how to under¬ 
stand the Z 2 invariant using the analogy of the global SU ( 2 ) anomaly 
[72], which reveals the relation between the second Ghern number and 
the mod 2 3d topological WZW term, as we now discuss. 

As discussed previously, see §3.4, the 2nd Ghern number C 2 for a time 
reversal invariant model is basically determined by the integral of the 
2nd Ghern class over If we decompose 3“^ into two hemispheres, 
then C 2 can be computed as the winding number over S^, 


( 77 ) 


C2 ^ 
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where g : ^ SU{2) is a transition fnnction on the overlap of two 

hemispheres (i.e. S^), for details see [56]. The right hand side is the 
winding nnmber of g aronnd the 3-sphere, i.e., the topological WZW 
term, and C 2 can be viewed as the obstrnction to the gange hxing 
problem. 

Taking the phase ambignity of a Kramers pair |\h) into considera¬ 
tion, when writing A = as explained in [56], the topological 

space where the transition fnnction really is dehned on is x S^. To 
accommodate the gange transformation r for the phase change, one 
needs an extra factor of SU{2), dne to an extra symmetry this action 
actnally lands in f/(l) C SU{2). It is natural, however to embed both 
these factors into 517(3), see below. Also, notice that the boundary 
d{S^ X D'^) = X and following Witten, the gauge held r can be 
extended to as follows. 


~g{x)=(^^^'^ .)eSU{3); r{9,p)=l G 5f/(3) 

^ \ -\/l - pe"*® / 

(p = 1 gives the 17(1) embedding for r.) The transition function is then 
given by 

gr ■■ X D‘^ ^ SU (3); (t, 9, p) ^ p)g{x)f{9, p) 


and the 5d topological WZW term T{gr) is dehned as usual. 


r(^.) 


—i 

240 ^ 




tr{gr ^dgrf 


which only takes values that are a multiple of tt, i.e., r(pr) G Z ■ vr. 
Actually, after integrating out the extra variables 9 and p, we recover 
the second Chern number again, 

(78) T{gr) = ^ J^ tr{g-^dgf = vrca 

We stress that the 5d WZW term is only well-dehned up to a mul¬ 
tiple of 27r because of the dependence of embeddings. In general, the 
embedding to 517(3) is not unique. Another possible way is to start 
with the mapping 

G : 5^ X 5^ ^ 517(2); {x, 9) ^ r-\9)g{x)r{9)g-\x) 

where r(9) is the diagonal matrix r{9) = clzap(e“*^/^, Since 

7ri(517(2)) = 0 and G(x, 0) = G(x, 27r) = 1, the map G : 5^ x 5^ —)■ 
SU (2) factors through 5^ in this case, analogously to §3.4. As a result. 
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G falls into two classes according to 7r4(S'17(2)) = Z 2 . Now extending 
G as before, we obtain a map G from into SU{3), 

G{x, 0, p) = gr{x, e, p)g;^{x, 0, p) 

It was proved in [34] that the corresponding WZW term is the same 
as before 

r(G) = ra) 

If we further choose a different 5-disc D'^ that bounds the same 
and then we consider the difference 

rD=(G) - rD«(G) = ^ j tr(G-HGf 

The right hand side is a multiple of 2?! because of 7 i^{SU{3)) = Z. 
This explains why we have to modulo 27i to eliminate the embedding 
ambiguity. Hence the 5d WZW term is well-dehned, i.e., independent 
of the embedding, modulo 27i, and by the relation (78) the 3d WZW 
term is then Z 2 -valued, 

(79) r53(5f) = 0,7r mod 27r 

In other words, the mod 2 second Chern number computes the topo¬ 
logical Z 2 invariant after embedding SU{2) x U{1) C SU{2) x SU{2) 
into SU{3). 

Summarizing the results: the phase ambiguity enlarges to S^x S^, 
which can be replaced by 5^ when a special boundary condition is 
satished [34], With the hxed target space SU{2), more room is made 
in the parameter space and the topological invariant originally assumed 
to lie in 7r3(S'f/(2)) ~ Z is now in 7 i4^{SU{2)) ~ Z 2 . In this language, a 
strong Z 2 topological insulator is characterized by a transition matrix 
(36) with odd winding number. 

4.3. Summary. It is shown in [47, 69] that for a 3d topological in¬ 
sulator, the Chern-Simons invariant and the Kane-Mele invariant are 
equivalent, i.e., v = u, since they are just the integral form and the dis¬ 
crete version of the mod 2 degree. Hence we call it the Z 2 invariant of 
time reversal invariant topological insulators and view it as a topolog¬ 
ical index. Another proof of the equivalence of the Chern-Simons and 
Kane-Mele invariants was given in [31] based on the group of reduced 
topological phases, which will be discussed in the next section. 

Earlier, we discussed a mod 2 index theorem, proved in [35], connect¬ 
ing the analytical index ind 2 {D) (50) and the topological index (38) 
for a Dirac operator with time reversal and chiral symmetries. View¬ 
ing this in the more general framework above, the reason why a mod 
2 index theorem appears in Z 2 topological insulators is the following. 
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In an odd dimensional time reversal invariant fermionic system, the 
parity anomaly pops up as a global anomaly, which is really difficult to 
compute in general [2] . By translating the problem into a gauge theory, 
the global parity anomaly is equivalent to a gauge anomaly, which can 
be dealt with locally. In our case, the parity of the spectral flow of edge 
states characterizes the global property of the material, and the gauge 
theoretic WZW term provides a practical way to compute the parity 
anomaly locally. It is well known that an index theorem can be used 
to compute a global quantity by a local formula, for example, the odd 
index theorem computes the spectral flow by the odd Chern charac¬ 
ter as discussed in §4.2.2. Therefore, the topological Z 2 invariant can 
be interpreted as the mod 2 version of the odd index theorem, which 
computes the global parity anomaly locally. 

Let us recap the relations between different variants of the Z 2 in¬ 
variant discussed above. Based on the non-abelian bosonization, the 
Z 2 invariant describes the parity anomaly of Majorana zero modes in 
a time reversal invariant topological insulator. On the one hand, the 
analytical Z 2 index is the parity of the spectral flow of edge states un¬ 
der an adiabatic evolution. The fermionic path integral of Majorana 
zero modes such as Dirac cones introduces the Pfaffian, while the path 
integral of the equivalent composite boson gives the square root of the 
determinant. After an adiabatic procedure, one can have a different 
way to bound two Majorana zero modes, so that the ratio of the effec¬ 
tive action, i.e., the sign of the Pfaffian, would change accordingly. The 
Kane-Mele invariant keeps track of the mod 2 version of the change of 
the signs of Pfaffians, that is, the parity of the spectral flow of chiral 
edge states through the adiabatic evolution. For inhnite dimensional 
Hilbert spaces, instead of the determinant (resp. Pfaffian) of matrices, 
we consider the determinant (resp. Pfaffian) line bundles of Dirac op¬ 
erators [29]. In the language of bundles, the Z 2 invariant is obtained 
by comparing the orientations of the Pfaffian and the determinant line 
bundles over the hxed points of the time reversal symmetry. If we 
further consider the holonomy of the Pfaffian or the determinant line 
bundle, the quotient form of the Kane-Mele invariant is replaced by 
a difference form, since the holonomy of the Pfaffian or determinant 
line bundle is an exponentiated eta invariant with a multiple of ivr [19]. 
If we model the edge states by Lagrangian sub-manifolds in the sym- 
plectic setting, then the spectral flow can be computed by the Maslov 
index, that is, the Lagrangian intersection number of the chiral edge 
states with the Fermi level. 
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On the other hand, the topological Z 2 index is the topological Wess- 
Zumino-Witten (WZW) term modulo two, since the bosonization pro¬ 
cess delivers the WZW term as the action functional of chiral edge 
currents [73]. The bulk-edge correspondence in this context is given by 
the Chern-Simons/WZW duality, and the WZW term is just the topo¬ 
logical index of odd Chern characters of gauge transformations. In our 
case, the time reversal symmetry induces a specihc gauge transform, so 
the Z 2 invariant is the mod 2 index of the odd Chern character of this 
special unitary in the odd K-group. Hence the mod 2 version of the odd 
index theorem computes the spectral flow of the edge Dirac Hamilton¬ 
ian, so the topological Z 2 invariant is interpreted as the Witten index 
Tr(—1)^. If we use the Feynman propagator instead of gauge trans¬ 
formations, the WZW term can be calculated based on the method of 
Green’s functions. The homotopy theory of Dirac Hamiltonians tells us 
that the 2nd Chern number, i.e., the integral of the 2nd Chern charac¬ 
ter is the fundamental invariant characterizing time reversal invariant 
models. In addition, the homotopy theory on the effective Hamiltoni¬ 
ans gives us an interesting explanation of the Z 2 invariant. According 
to the dimensional ladder in anomalies [56], the 3d parity anomaly is 
the descendant of the 4d chiral anomaly, so the WZW term (its mod 2 
version is the parity anomaly) is naturally connected to the 2nd Chern 
character (which computes the chiral anomaly). In other words, the 
topological Z 2 invariant can be derived from the 2nd Chern number 
from the relation between anomalies. Finally an analogy of the global 
SU{2) anomaly of the WZW model explains that it is naturally Z 2 - 
valued, so we identify the topological Z 2 index as the mod 2 WZW 
term, and it is called the parity anomaly instead of the SU{2) anomaly 
in condensed matter physics. 

5. K-theoretic classification 

In this section, we will discuss the K-theoretic classihcation of non¬ 
interacting (or weak interacting) Dirac Hamiltonians of topological in¬ 
sulators. More precisely, the Z 2 -valued topological invariants of topo¬ 
logical insulators and Bogoliubov-de Gennes (BdC) superconductors 
fit into a periodic table resembling the Bott periodicity as in topo¬ 
logical K-theory. First, a Clifford algebra classification for these was 
established in [62] and the K-theory classification followed as the Clif¬ 
ford extension problem [51]. A hner classihcation of topological phases 
based on KR-theory and twisted equivariant K-theory was proposed in 
[31]. For time reversal invariant topological insulators, KH-theory was 
also applied to the classihcation problem in [26]. 
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KR-theory, introduced by Atiyah [4], is the K-theory of Real vector 
bundles over an involutive space. By contrast, KH-theory is the K- 
theory of Quaternionic vector bundles over an involutive space. They 
are isomorphic to each other by a Fourier transform. A detailed treat¬ 
ment of twisted equivariant K-theory can be found for example in [30]. 

5.1. Tenfold way. Generalizing the standard Wigner-Dyson threefold 
way, namely, symmetries in quantum mechanics described by unitary, 
orthogonal and symplectic groups, Altland and Zirnbauer [1] further 
proposed a tenfold way in random matrix theory, which implies that 
free fermionic systems can be classihed by ten symmetry classes. In¬ 
deed, there are three types of discrete (pseudo)symmetries in topo¬ 
logical insulators and BdG Hamiltonians: the time reversal symmetry 
T, the particle-hole symmetry V and the chiral symmetry C, and the 
combinations of these three symmetries give ten classes in total. 

Based on a time reversal invariant Dirac Hamiltonian H, which 
means that = H, and additionally = ±1 depending on the 

spin being integer or half-integer, the time reversal symmetry (TRS) 


induces three classes 



f +1 

if 

TH{k)T-^ = H{-k), = +1 

(80) TRS = { -1 

if 

rH{k)r-^ = H{-k), = -i 

1 0 

if 

rH{k)r-^ ^ H{-k) 


Similarly, the particle hole symmetry (PHS) also gives three classes, 

r +1 if VH(k)V-^ = -H(k), = +1 

(81) PHS = < -1 if VH{k)V~^ = = _i 

[ 0 if VH{k)V-^ 7^ -H{k) 

The chiral symmetry can be dehned by the product C = T ■ V, some¬ 
times also referred to as the sublattice symmetry. Since T and V are 
anti-unitary, C is a unitary operator. If both T and V are non-zero, 
then the chiral symmetry is present, i.e., C = 1. On the other hand, 
if both T and V are zero, then C is allowed to be either 0 (type A or 
unitary class, according to the classihcation of [1], see e.g. [62]) or 1 
(type AHI or chiral unitary class). In sum, there are 3 x 3 -|- 1 = 10 
symmetry classes. In particular, the half-spin Hamiltonian with time 
reversal symmetry falls into type AH or symplectic class, which is the 
case we are mostly interested in, for more details about the 10-fold way 
and other symmetry classes see [62]. 

In [51], the above particle-hole symmetry V is replaced by a U{1)- 
symmetry Q representing the charge or particle number. Roughly, 
the parity V is related to the charge hj V = (—1)®. So sometimes 


50 RALPH M. KAUFMANN, DAN LI, AND BIRGIT WEHEFRITZ-KAUFMANN 

topological insulators are also referred to as U{1) and time reversal Z 2 
symmetry protected phases. 

5.2. Nonlinear cx-model. One classification scheme is given by the 
homotopy groups of non-linear sigma models or target spaces of the 
Brillouin torus [62]. Assume that there are N occupied bulk bands 
and denote the occupied bands by |k) := Further as¬ 

sume that there are also N unoccupied bands and that there exists a 
band gap between the conduction (occupied) and valence (unoccupied) 
bands, and the Fermi level Ep is assumed to lie in the band gap. If 
there is a physical material the Fermi level is determined by it, purely 
mathematically, without a given material, one can choose Ep to be 
any energy in the band gap. Dehne the Fermi projection in ket-bra 
notations as 

N 

(82) Pf := 5]] |n,k)(n,k| = |k)(k| 

71=1 

In addition, the Fermi level is always assumed to be zero, i.e., Ep = 0, 
so that the occupied bands have negative eigenvalues. From now on 
we denote the Fermi projection by P = Pp for simplicity. 

For non-degenerate band structures, the spectral flattening trick can 
always be applied, which replaces the Pth band Pj(k) (a continuous 
function over the Brillouin torus) by a constant function Pi(k) = Pj. 
Dehne the sign of the Dirac Hamiltonian by P(k) := P(k)|P(k)|“^, by 
which the band structure is reduced to a battened two-band model, or 
in terms of the Fermi projection the battened Hamiltonian P is written 
as 

P(k) = 1 - 2P 

We use the notation P to indicate that it looks like the grading operator 
of a Fredholm module in K-homology satisfying 

p2 = F* = E 

For the unitary class (type A), the battened Hamiltonian P(k) debnes 
a map from the Brillouin torus to the Grassmannian, aka. a non-linear 
sigma model: 

(83) P : T'^ ^ Gr{N, C^^) = U{2N)/{U{N) x U{N)) 

since the occupied and empty bands each have N levels by assumption 
and we don’t distinguish energy levels among occupied or empty bands. 
Similarly, for the symplectic class (type AH), P debnes another non¬ 
linear sigma model, 

(84) P : ^ Gr{N, M^^) = 0{2N)/{0{N) x 0{N)) 
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Without loss of generality, we assume that the empty bands have 
inhnitely many levels. In other words, the 2iV-dimensional total vector 
space is replaced by the inhnite dimensional Hilbert space, then the 
classifying space BU{N) (resp. BO{N)) plays the role of the Grass- 
mannian Gr{N, C‘^^) (resp. Gr{N, If N is big enough, the lower 

dimensional homotopy groups stabilize. Accordingly one can regard the 
target spaces the direct limits, 

BU = \i^BU{N), BO = \\^BO{N) 

N N 

which capture the stable (under increasing N) features, which is what 
stable homotopy theory is concerned with. 

For a compact Hausdorff space X, the topological K-theory K{X) is 
dehned to be the Grothendieck group of stable isomorphism classes of 
hnite dimensional vector bundles over X. The complex or real K-group 
can be equivalently dehned by homotopy classes of maps from X to the 
classifying spaces, 

K{X) = Kc{X) ~ [X,BU X Z], KO{X) = K^{X) ~ [X,BO x Z] 

In stable homotopy theory, we say the spectrum of the complex K- 
theory is BU x Z, which has a period 2, i.e., 0?BU BU x Z, where 
G is the loop space operation. Similarly, the spectrum of the real K- 
theory is BO x Z, which has a period 8, i.e., Q^BO ~ Z x BO. Such 
periodic phenomena, are known as the Bott periodicity. The hrst case 
was discovered by Bott in the study of homotopy groups of spheres 
using Morse theory. For unitary or orthogonal groups, Bott periodicity 
is spelled out as 

TTk+iiBU) = TTkiU) ~ 1Tk+2iU), TTk+iiBO) = 71^(0) ~ TTk+siO) 


Besides in topological K-theory and homotopy groups of classical groups, 
Bott periodicity is also present in Glifford algebras and their repre¬ 
sentations. Of course, Glifford modules and topological K-theory are 
intimately related in index theory. 

Going back to band theory, in the language of K-theory via stabi¬ 
lizing (embedding into the direct limit) the non-linear sigma model 
induces an element in the complex or real K-group, the homotopy class 
of the map T'' A Gr(N) BU(N) -> BU. 


(85) 


rp, f iF(T'^) for type A 

^ ^ KOiT'^) for type All 


Furthermore, if the map from the Brillouin torus can transformed into 
a map F from the sphere either by a qnotient map or by extending 
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and quotienting, as described in the previous section, then 


( 86 ) 


f TTniBU) = TTn-iiU) for type A 

[ TTniBO) = TTn-iiO) for type All 


For instance, consider the edge Hamiltonian in the quantum spin Hall 
effect, i.e., n = 1, then the 2 d Z 2 invariant corresponds to 


7ri(HO)=7ro(0)~Z2 


5.3. Clifford modules. Following [51], we briefly go over the classih- 
cation approach via the Clifford extension problem in this subsection. 
For this one regards free fermions in the basis of operators Cn which 
satisfy the Clifford relations C/Cm + CmQ = 25/^m- 

Then the B(l) symmetry Q and conjugation by the time reversal 
symmetry T have matrix representations on the Cj 

(87) Q = foo ® o-y, T = —iay 0 

where (Tq = and cTj are the Pauli matrices. 

Q'=T2 = -1, {Q,T} = 0 

A Dirac Hamiltonian H = Aj^CjCk represented by a real skew- 
symmetric 2 Nx 2 N matrix A has U{1) symmetry if and only if [Q, A] = 
0. Similarly, H being time reversal invariant is equivalent to {T,A} = 
0. To simplify things, one can rescale A, e.g. by regarding it as a 
skew-symmetric quadratic form or iA as a Hermitian form, to have 
Eigenvalues i,—i. Denote the result by A, then = — 1 

Finding a compatible matrix A reduces to the problem of Ending an 
extension of a Clifiord algebra with n generators by one more generator 
A. Here n = 0 for no symmetry, n = 1 for T only and n = 2 for T 
and Q. For example for n = 2 : ci = T and 62 = QT and we want to 
add one more generator 63 = A for the Dirac Hamiltonian with both 
time-reversal symmetry and U{1) symmetry. In general, the Clifiord 
extension problem is given by considering the extensions 

i '■ C/o,n(I^) ^ C/o,n+l(IR) 

Here Clp^q is the real Clifiord algebra for the quadratic form with sig¬ 
nature (+,..., -F, -,..., -). 

-V-^ -V-^ 

p g 

Bott periodicity in Clifiord algebra manifests itself as Morita equiv¬ 
alence, 

Cln+2{C) ~ M2{CL{C)) C'fo(C), 

Clp,q+s(^) — Clp^8^q{M.) ~ Mie(C/p_g(M)) Clp^q{M.) 
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where Mn{R) is the n x n matrix of some commutative ring R and 
~M stands for the Morita equivalence in ring theory. Recall that two 
unital rings are Morita equivalent if they have equivalent categories of 
left modules. 

In order to classify the Dirac Hamiltonians with or without the above 
mentioned Z 2 symmetries, we actually deal with the representations 
of the above Clifford algebras, i.e. Clifford modules, since the are 
represented by matrices. Denote by M (C'/o,n(I^)) fhe free abelian group 
generated by irreducible Z 2 -graded C'/o,n(l^)-iiiodules. There exists a 
restriction map induced by z. 

Each Clifford extension z defines the cokernel of i*, 

(88) := cokerit*) = M(C/o,n(M))/z*M(C'/o,n+i(M)) 

For the complex case, we dehne M‘^(C/o,n(l^)) as the Clifford modules 
of the complexihcation Cln{C) = Clo^niM) C)r C, can be obtained 
similarly. It is An (or A^) that encodes the information about non¬ 
trivial extensions of free Dirac Hamiltonians of different dimensions in 
the periodic table. 

Furthermore, the Clifford modules and the cokernels inherit the Bott 
periodicity: 

Ak+8 ^ Ak, Hfc+2 - 

In fact, the cokernels have the same form as the homotopy groups of 
the classifying spaces [5], 

An-TTniBO), ~ 7r„(Rf/) 

For instance, if we extend C'/o,i(K) generated by T in the quantum spin 
Hall effect, 

z : C'Vi(M) ~ C ^ Clo, 2 {^) ~ M 

then Ai ~ Z 2 gives the 2d Z 2 invariant. One step further, if we extend 
C/ 0,2 generated by T and QT for 3d topological insulators, 

i : C/ 0,2 (M) ~ e ^ C/ 0,3 (M) ^ H © e 

then A 2 — Z 2 shows that there exists another Z 2 invariant in 3d. 

It is easy to promote vector spaces to vector bundles and get the 
K-theory of Clifford bundles, that is, to pass from the Clifford algebra 
classification to the K-theory classification of topological insulators. 
For any pseudo-Riemannian vector bundle E ^ X over a compact 
Hausdorff space X, dehne an associated Clifford bundle Cl{E) —)■ X as 
usual, and a representation gives a Z 2 -graded vector bundle M{E) —)■ 
X. The Clifford extension problem arises naturally when we add a 
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one-dimensional trivial bundle Ti —)• X by the Whitney sum, and the 
cokernel bundle is dehned as 

A{E) := cokerii* : M{E © Ti) ^ M{E)) 

By the Clifford grading M{E) = and its Euler class is 

defined as the formal difference 

Xe{M) := [M°] - [M^] 

which induces a map from the cokernel bundle to the reduced KO- 
group, ^ 

XE : A{E) ^ KO{B{E), S{E)) = KO{E+) 
where the right hand side is the Thom construction, B{E) and S{E) are 
the ball and sphere bundles of E, E~^ is the one-point compactification. 
In particular, if the vector bundle E is of rank n, and assume the base 
space is contractible X ~ pt, then A{E) is reduced to An and we get 
the isomorphism 

A, ~ KO{S^) = TTniBO) 

More details can be found in [5]. 

On the level of Clifford modules, the Clifford extension problem is de¬ 
scribed by the triple {E, E; a), where E, E are elements in M{CIq ,rt+l(l^)) 
and cr is a linear orthogonal map identifying E\ciq^„ = -E|c/o,„ as re¬ 
stricted representations of C/o,n(lR), which has a direct generalization 
into Clifford bundles. In general, for a pair of spaces {X, Y), given two 
vector bundles E,E on X and a bundle map a : E ^ E such that a 
is an isomorphism E\y — E\y when restricted to Y, the Atiyah-Bott- 
Shapiro construction [5] defines an element [{E,E]a)] in the relative 
K-group K{X,Y) := K{X/Y). In [51], Kitaev applied the Atiyah- 
Bott-Shapiro construction to the Clifford extension problem on Clifford 
bundles and proposed the K-theory classification scheme for topological 
insulators. 

5.4. Twisted equivariant K-theory. An even more general frame¬ 
work has been set up in [31] for studying topological phases of quantum 
systems with additional symmetries by twisted equivariant K-theory. 
For a symmetry group G, a homomorphism 0:G—)-{±l}is defined 
to keep track of the unitary or anti-unitary symmetry, and a group 
extension r, 

1 ^ T ^ G" ^ G ^ 1 

is introduced to take care of the phase ambiguity in quantum mechan¬ 
ics. As a topological invariant, the group of reduced topological phases 
RTP{G,(f),T) is defined to be the abelian group of generalized quan¬ 
tum symmetries {G,(j),T). In this subsection, we will see how to apply 
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twisted equivariant K-theory to classify topological phases closely fol¬ 
lowing [31]. 

For time reversal invariant topological insulators, the reduced topo¬ 
logical phase group can be computed by the shifted KR-group of the 
Brillouin torus, 

(89) RTP{G, 0, r) ~ KR-^T’^) 

Similarly, if the parity is reversed (i.e., with particle-hole symmetry V 
and = — 1), it can be computed by the equivariant K-group, 

(90) RTP(G,0,r) ~iFz,(T'') 

Finally, if = —1 and = —1, then the twisted equivariant K- 
groups are applied, 

(91) RTP(G, 4>, t) ~ ~ A-Oi*(T'') 

For the torus T'^, the KR-groups can be computed by the stable 
splitting of the torus into wedge products of spheres and the twisted 
Thom isomorphism. Using the reduced KR-theory dehned as usual 
KR(X) ~ ^(X) © Z, we have 

- ^ 2 , 

Let us look closely at two natural maps appeared in the computation 

of KR (T^). One is the projection map, i.e., three projections Pij : 
—)■ T^, which induce three embeddings, 

i.e. Z2-A (Za)^^ 

The other one is the collapse map g : —)■ 5^, which also induces a 

map between the KR-groups, 

q* : ^ 

Further, it was shown in [31] that its image gives the fourth Z 2 com¬ 
ponent in KR (T^), 

/m(g*) ~Z2 -A (Z2)^" 


which corresponds to the Z 2 invariant for strong topological insulators. 
First note that the stable splitting of is given by 

^stable 


'-4 


The idea behind the proof is that a canonical element in KR (S^) is 
obtained by constructing a principal S'p(l)-bundle E on S^, then pull 

it back to get an element q*E G KR (T^). Taking the time reversal 
symmetry into account, the image q*E is then identihed with the fourth 
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'-4 , 


component Z 2 G KR (S^) from the stable splitting. In the 2d case 
the situation is similar and actually one can prove which component is 
used by looking carefully at the construction of §3.4. 

The group of reduced topological phases is important in that it pro¬ 
vides a universal topological invariant. Recall that u is the Kane-Mele 
invariant dehned by the change of the signs of Pfaffians, and v is the 
Chern-Simons invariant dehned by the mod 2 WZW term. For 3d 
topological insulators, Freed and Moore [31] gave another proof of the 
equivalence u = v hy showing that both u and v factor through the 
abelian group RTP(G, 0, r) = KR~'^{T^), furthermore both vanish on 
Im{p*j) and equal the identity on Im{q*). 


6. Further Discussions 

6.1. Bulk-edge correspondence revisited. In this subsection, let 
us go over the bulk-edge correspondence in topological insulators again 
on different levels. We will relate it to a generalized index theorem 
for further investigation. In general, the bulk-edge correspondence is 
a realization of the holographic principle in condensed matter physics. 
The philosophy of the holographic principle is that the boundary held 
theory determines the geometry in the bulk, and conversely the bulk ge¬ 
ometry determines the boundary theory up to some equivalence. Even 
though the general holographic principle still remains a conjecture, 
the Chern-Simons/Wess-Zumino-Witten (CS/WZW) correspondence 
is well understood as an actual theorem, which actually gives the bulk- 
edge correspondence for 3d topological insulators on the action func¬ 
tional level. 

On the level of effective Hamiltonian, the Id edge Hamiltonian is 
obtained from the 2d bulk Hamiltonian by a partial Fourier transform 
along the direction with translational invariance in §2.4. Such corre¬ 
spondence realized by the partial Fourier transform holds for any di¬ 
mension, so there exists a one-to-one correspondence between the bulk 
and edge Hamiltonians. 

For 3d topological insulators, the Z 2 invariant is the mod 2 WZW 
term, which is derived from the Chern-Simons theory in §4.2.1. In 
other words, the bulk-edge correspondence on the level of action func¬ 
tional is given by the duality between the bulk Chern-Simons action 
and the boundary topological WZW term. If we consider quantum held 
theories, then the bulk is described by a topological quantum held the¬ 
ory (TQFT), and the boundary is characterized by a rational conformal 
held theory (RCFT). So the bulk-edge correspondence on the level of 
held theory is given by CS3/WZW2 as a correspondence between 3d 
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TQFT and 2d RCFT. On the other hand, based on the gauge group 
G such as SU{2), there exists a one-to-one correspondence between 
the WZW models and Chern-Simons theories, both characterized by 
their level fc G Z. So on the level of group cohomology, the bulk-edge 
correspondence is given by the Dijkgraaf-Witten theory [27]. If the 
Dijkgraaf-Witten transgression map is r : —)■ H^{G,'L), 

then T is the correspondence between the Chern-Simons theory and 
WZW model associated with G, for the geometric construction on bun¬ 
dle gerbes, see [22], 

The generalized index theorem in §2.4 can be viewed as a bulk-edge 
correspondence on the level of index theory. With the presence of a 
domain wall, two quantum Hall states characterized by distinct hrst 
Chern numbers are separated by a boundary created by the domain 
wall. The variation of the hrst Chern numbers across the domain wall 
is identihed with the spectral how of the edge Hamiltonian along the 
boundary. In this model, the index map of the extended Hamiltonian 
plays the role of a correspondence connecting the topological invariants 
from the bulk and the analytic behavior of the edge states. 

For 3d Z 2 topological insulators, we have a similar picture as that in 
the above quantum Hall model. With the presence of the time reversal 
symmetry, two 3d bulk systems characterized by distinct Chern-Simons 
theories are separated by a boundary consisting of the hxed points of 
the time reversal symmetry. If the time reversal symmetry is encoded 
into a specihc gauge transformation, then the variation of the bulk 
Chern-Simons actions is exactly the topological WZW term of this 
specihc gauge transformation. By the global SU (2)-anomaly, the WZW 
term is indeed Z 2 -valued. On the other side, the Pfafhan formalism 
of the Kane-Mele invariant counts the parity of the spectral how of 
the chiral edge state along the boundary, i.e., the hxed points. Now 
the correspondence, i.e., the index map, is given by the Witten index 
Tr(—1)^ since it counts the parity of Majorana zero modes such as 
Dirac cones. Hence the story of the topological Z 2 invariant can be 
interpreted as a bulk-edge correspondence realized by the Witten index. 

If we model the bulk and boundary spaces by G*-algebras, then 
the bulk-edge correspondence can be realized by Kasparov module in 
bivariant K-theory. So on the level of Kasparov’s analytic bivariant 
K-theory, the bulk-edge correspondence is the Kasparov product with 
a carefully selected Kasparov module according to the interpretation 
of index maps in KK-theory. For quantum Hall ehect, such a bulk- 
edge correspondence is explicitly constructed in [20] , where the authors 
represent the bulk conductance as a Kasparov product of the boundary 
K-homology with a specihc extension class (i.e. KK^-c\ass) as the 
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correspondence. A similar model for time reversal invariant topological 
insulators is expected to exist, and the bulk-edge correspondence is 
expected to be realized by Kasparov product in real KK-theory. 


6.2. Noncommutative geometry. There are different sources of non¬ 
commutativity in condensed matter physics, which call for generaliza¬ 
tions of classical models in the noncommutative world. For instance, 
the introduction of magnetic helds into a lattice model breaks the trans¬ 
lational symmetry, that is, the effective Hamiltonian does not commute 
with the translation operators any longer, which is the key assumption 
for the classical Bloch theory of crystal wavefunctions. Another source 
of noncommutativity is created by nontrivial quasi-particle statistics 
such as non-abelian anyons in fractional quantum Hall systems. In 
addition, systems with impurities or defects such as disordered topo¬ 
logical insulators are of realistic importance, which is the main source 
of noncommutative geometry models up to date. 

Noncommutative geometry [23] has already achieved some success 
in modeling the integer quantum Hall effect and topological insulators 
with disorders. The work by Bellissard, van Elst and Schulz-Baldes [16] 
set up the framework for studying topological invariants of a disordered 
system by noncommutative geometry. The disordered integer quantum 
Hall effect can be characterized by the noncommutative hrst Chern 
number, which was obtained by combining the Connes-Chern character 
with the Kubo-Chern formula [16]. The noncommutative 2nd Chern 
number is a natural generalization for time reversal invariant models, 
and the general noncommutative n-th Chern number was also studied 
in [58]. 

A survey on the noncommutative approach to disordered topological 
insulators was written by Prodan [57]. The noncommutative odd Chern 
character was discussed in a recent paper [60] by Prodan and Schulz- 
Baldes in the same framework as in [16]. A noncommutative formula for 
the isotropic magneto-electric response of disordered insulators under 
magnetic helds was also proposed in [53]. It would be interesting to 
see more noncommutative generalizations of both the analytical and 
topological Z 2 index of topological insulators in the future. 
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